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1. Introduction 



Tunneling phenomena are unique to quantum systems and play important roles in various 
fields. To name a few of them, B+L violation of the standard electroweak model, inflation in 
the early universe, nuclear fission and fusion reactions, spontaneous nucleation in He 3 /He 4 
systems, and chemical reactions are the area of immediate influence. Although it can be 
readily understood in Schrodinger formalism in simple quantum mechanical systems, once 
the degree of freedom increases the calculation becomes a challenge. This situation is quite 
serious in field theories, which have infinite degree of freedom. Since the tunneling phe- 
nomena is essentially non-perturbative one, extracting relevant degree of freedom and doing 
reliable calculation requires new formalism. 

The imaginary time path-integral method has been known to be effective in dealing with 
such a situation for some time. This is mainly because of the existence of the classical con- 
figurations, instantons and bounces, around which the path-integral may be evaluated. In 
theories with degenerate classical vacua, one often has instanton solutions, which connects 
the different (perturbative) vacua. When the coupling is small enough, the dilute-gas approx- 
imation for multi-instanton configurations yields the energy splitting of the ground states by 
the tunneling. Bounces are configurations that starts from the ground state and "bounces 
back" at the escape point and then comes back to the original ground state. A relevant ana- 
lytical continuation of the Gaussian integral around bounce solution yields the decay rate of 
the unstable states. The power of the instanton method and the bounce method described 
so far are actually quite limited. Here we will point out three of the problems encountered 
by these methods. 

The most obvious one is the fact that those calculations are limited only to the ground 
state. In quantum mechanical context, if one wishes to obtain the energy splitting of the 
excited states in the symmetric double-well, or if one wants the decay rate from higher 
excited states, one realizes the current methods are not enough. In field theoretical context, 
this difficulty became most apparent and serious in the late 80's, when the baryon number 
and lepton number violation processes in the electroweak standard model caught attention 
of some of the particle physicists, in relation of the anticipated PP collider in TeV region 
at that time. The tunneling process induced by the colliding quarks results in the B+L 
violation through the anomaly of the model. Since the potential height is of order M\y /a ~ 
few TeV, this process is highly suppressed in the low energy region. However when the initial 

*' In some literature, the word "instanton" is used for any the imaginary time classical configurations. 
This, however, is avoided in this paper for the sake of clarity; "instanton" and "bounce" are distinguished 
as will be seen in the rest of this paper. 
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energy of order TeV is available, this process may occur at observable rate. However, it was 
soon noticed that the usual instanton calculation is not only invalid, but yields inconsistent 
result at higher energies, violating unitarity. Thus the necessity of the further development 
of the formalism became apparent. 

Another difficulty is a subtle one. There are theories in which tunneling is physically 
possible, but no classical solution exists. This happens for a specific potential in the scaler 
field theories and in the gauge-Higgs system. An elementary scaling argument shows that 
the action of any finite action configuration can be lowered by shrinking it. In other words, 
any configuration that is the stationary point of the action functional is point-like, which 
denies any hope of naive application of the instanton formalism. 

The third problem is actually a paradox. In the quantum mechanics with asymmetric 
double-well potential, the ordinary bounce method yields a decay rate. This is due to the fact 
that the Gaussian integration of the fluctuations is concerned only with the neighborhood of 
the bounce solution. On the other hand, there is no decay, since Hamiltonian is hermitian; 
the energy spectrum is certainly affected by the tunneling, but does remain real and discrete. 
Thus this instability is a "fake" . 

Motivated by these problems at various stages, we have developed several formalisms 
in the recent years. The aim of the paper then is to present these developments in a self- 
contained manner, together with several new results. We believe that the power of these 
formalisms is not limited to the original contexts and thus hope to make these widely available 
to researchers in other fields as well. 

The organization of this paper is as follows. In section 2, we briefly review the B+L 
violation problem in the standard model, which was the starting point of these works. The 
use of the dilute instanton-gas approximation, its infamous "unitarity problem" , its solution 
by the "interacting instanton method" , the need for further developments are explained in a 
concise manner. The ultimate method in this line is "the proper valley method", described 
in section 3. The presentation here is in a general model, not limited to any particular 
quantum mechanical or quantum field theoretical model. Its advantage over the pre-existing 
methods, the streamline method and the constrained instanton method are explained in 
detail. A numerical analysis that demonstrates the power of the proper valley method is 
also presented. 

The application of the proper valley method to various models are presented in section 4. 
The asymmetric double-well potential, the "Fake instability" problem is treated and solved 
in the proper valley method in subsection 4.1. False vacuum problem in the scalar field 
theory is analyzed in subsection 4.2. Unstable scalar field theory that does not possess any 
finite-size instanton is studied in subsection 4.3, where we also compare our results with 
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the one from the constrained instanton method. Similar analysis is done for the gauge- 
Higgs system in subsection 4.4. Subsection 4.5 gives the application of the valley instanton 
configuration obtained in subsection 4.4 to the actual calculation of the cross section in the 
standard model. 

In section 5 we describe a completely different approach, the complex time method. The 
big advantage of this formalisms over the imaginary time formalism is that the configuration 
goes through not only the forbidden region (just like instanton/bounce), but also the allowed 
region. This implies the solution of the fake instability problem, and calculation of the excited 
states, among others. However, the difficulty lies in the fact that the exact formulation is not 
clear at all. It is easy to construct complex-time configuration by patching up the instanton 
and the real-time solutions and then use those for evaluation by hand- waving argument. 
However, just how exactly these configuration should contribute to the path- integral remains 
a mystery, even in a simple one-dimensional quantum mechanics. What we present in this 
section is an analysis of the WKB calculation to show that there is at least one formalism 
that works in the one-dimensional quantum mechanical system. 

In section 6, we discuss the problems of the perturbation theories that arise due to the 
existence of the tunneling phenomena. In models whose the perturbative series is not a 
Borel-summable series, we expect that it is due to the fact that physical observables are 
well defined only after the non-perturbative terms are added. This interplay between the 
perturbative effects at large orders and the non-perturbative effects has been a subject of 
the investigation for some time. In this section we give some of the new insights obtained 
recently. In subsection 6.1 we give a Borel analysis of this phenomenon using the idea of 
valley. And then in subsection 6.2 we present the theory of the asympton, which is the 
configuration that fictitiously dominates the perturbative functional at large orders. This 
provides us with an insight as to the original of the non Borel-summability and the possible 
direction for the converging calculations. 



§2. Interacting instantons in the anomalous B and L violation 

The anomalous baryon number (B) and lepton number (L) violation is one of the most 
fascinating processes predicted in the electroweak standard model. The prediction is based 
on the chiral anomaly and topological transition. In the electroweak model the apparent 
conservation of baryon number current jg and the lepton number current j£ is, in fact, no 
longer valid due to the anomaly. Their divergence is deduced from the representation of the 
quarks or leptons under SUl(2) and Uy(l) group. Since we have, in each generation, three 
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SU(2) doublets that has Y = 1/3 for the left-handed quarks and six singlets with Y = 4/3 
or —2/3 for the right-handed ones, the divergence of jg is 

Nt 

d,f B = ^- 2 [9 2 w; y W a ^ - g' 2 B, u B^\, (2-1) 

where 

W; v = W - d v W; - ge abc WlW c ul (2-2) 

B^v = d^B v — dyB^ (2-3) 

are the field strength of the SU(2) W® and U(l) B^ gauge fields, respectively; g and g' are 
their coupling constants; W afU/ = (l/2)e flvXa W% al and the iVf is the number of the generation. 
Interestingly the divergence of j£ is exactly the same. Therefore the (B— L) current stays 
conserved. 

The topological transition resides in the SU(2) sector in the model. Consider a trajectory 
in the finite energy configuration space made of WHx) and the doublet Higgs <P(x) that starts 
from and ends up with the vacuum. We parametrize this as 

W;(0, x) = W%(T, x) = 0, <Z>(0, x) = $(T, x) = 

using a parameter t G [0,T], where v is the vacuum expectation value. All the trajectories 
of this type are divided into topological classes that are labeled by the winding number of 
the homotopy n 3 (S 3 ). The integral along the trajectory, 

2 

Nw ^3hJ dt d ' X W l^\ (2-4) 

gives its winding number. Combined with Eq. ( |2-1| ), this relation indicates that the transition 
along the trajectory causes the B and L violation by N{N W units. 

The height of the potential barrier that the topological transition with iV w = ±1 need 
to go over is of the order of (47r/g 2 )mw This is the static energy E sp of the sphaleron, a 
static but unstable solution of the equations of motion. 0) This large energy causes a strong 
suppression if the transition takes place through quantum tunneling starting from the low- 
lying states near the vacuum. Actually the WKB suppression factor in the amplitude is 
estimated!) to be e - 8n2/9 ' 2 ~ e" 180 . 

A naive expectation, however, leads us to the interesting possibility that this strong 
suppression may be weakened if the energy of the transition becomes comparable to E sp . 
Along this line, two issues have been discussed in literature; the effect of the B and L 
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violation on the baryogenesis of the Universe Ef and the direct detection of the violation at 
a high energy collision. & 

Only reliable estimates of the transition amplitudes in high energy processes can answer 
the question about these issues. This is not, however, a simple task. The difficulty lies in 
the fact that it is non-perturbative: the process 

q + q-*7q + 3l, (2-5) 

or, more generally, the one accompanied by the simultaneous production of n w W bosons 
and 7Zh Higgs bosons, 

q + q -> 7q + 3l + n w W + n u 0, (2-6) 

is the specific example of the transition, but we cannot write down the relevant diagrams for 
this process by using the Feynman rules in the electroweak standard model. 

't Hooft is the first person to have estimated the amplitude. & His tool was the imaginary 
time formalism of the path-integral; it can incorporate the topological transition as the 
background in the path-integral. Ringwald and Espinosa pushed the idea further; they 
applied the path-integral to the evaluation of Green functions appropriate for the scattering 
process ( |2-6| ).li'll§ Aoyama and Kikuchi improved their analysis by taking into account multi- 
transitions; EP this can resolve the unitarity problem pointed out in their single transition 
calculation. The following of this section is a review of this non-perturbative calculation 
according to the line mentioned above. We are still at the middle of the way to the final 
answer and a lot of related works have been published. We are not going to look back all of 
these works here. (For the other related works, see, for example, Ref. ^).) 

Let us start with Ringwald and Espinosa calculation. For simplicity, we take the SU(2)- 
Yang-Mills-Higgs-model with a single generation of massless quarks and leptons. All the 
essence of the calculation can be seen in this simple version. We intend to evaluate the B 
and L violating Green functions of the Heisenberg operator, 

G(x\, £4; yi, y nw \ z ii •••) z n H ) 

_ (yac\q 1 (x 1 )q 2 (x2)q3(x3)l(x4)W(y 1 )..W(y nw )(f)(z 1 )..4(z nH )\veLc) 

(vac I vac) 

where qi(x), (feOz), and q${x) are SU(2) doublet quarks distinguished by their color, l(x) the 
doublet lepton, <f>(x) the physical Higgs boson. We have suppressed the Lorentz and the 
SU(2) indices in Eq. (2-7). We use the imaginary time formalism for the calculation, i.e. 
the time evolution operator is e~ rH (H is the Hamiltonian) instead of e~ ltH . This Green 
function can be evaluated also in the path-integral 

G(xi, 24; y nw 'i Zi, z nH ) 
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= 4/ II Wdf\[dW][d#\exp{-S\W,$,q,l\} 
f=qi,q2,qs,i 

xqi(x 1 )q 2 (x 2 )q 3 (x 3 )l(x i )W(y 1 )...W(y nw )^(z 1 )...^)(z nH ), (2-8) 

Z = J J] [dfd/][dW][^]exp{-^^g,Z]}, (2-9) 

/=?l,92,93,i 

where S 1 is the sum of the terms Sw> S^, and Sf , evaluated respectively from the Lagrangians, 

^ = \(W« V )\ (2-10) 

C^ = {D^ 2 + \(\<P\ 2 -v 2 ) } (2-11) 

A = QcD^a^qc + ID^a^l. (2-12) 

c 

In these equations, is the covariant derivative; 



D, = d, + ig[-)W^ (2-13) 



and cr M = (a, — i) is the two component Pauli matrices for the left-handed spinorsEll. Here- 
after, the Greek indices are understood as being summed up in the Euclidean metric. 

The usual perturbative calculation of the function (2-8) around the vacuum is trivially 
zero. We instead evaluate it in the background of the topological transition. The anomaly 
equation (|2-1|) indicates it should have |iVw| = 1. In the standard model, we do not have a 
stationary configuration of the action around which we construct an systematic expansion 
of the Green function in powers of g: any finite size configuration tends to shrink toward 
zero size. To circumvent this, we use the constrained instanton method. We first restrict 
the functional space by imposing the constraint on the size p and evaluate the path-integral 
over all fluctuations around the stationary configuration in this restricted functional space. 
At the end of the calculation, we integrate over p to recover the whole functional space. 
The constrained instanton is the name for the topological transition stationary in the re- 
stricted functional space. Hereafter in this section we simply call it instanton. The bosonic 
component of the instanton in the singular unitary gauge has the following form 

Wi^-^Alnfl + ^j, (2-14) 



*> We use the relation t = —it of the imaginary time r and the real time t when we need to go back to 
Minkowski expression. Using the other one t — it as in Refs. is just the matter of convention. 
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in the central region x <C m^, rrig 1 , and 

W in ; ~ ^Lf-ria^dvG^x), (2-16) 

(P Jn ~(l-27r 2 ^G mH (x)) ^ ° j, (2-17) 

in the asymptotic region x p, where we have used G m for the four dimensional massive 
propagator, 

(-d 2 + m 2 )G m (x) = 5\x), (2-18) 

and the 't Hooft symbol r] afMU , 

' 5 av if /I = 4; 



if i/ = 4; (2-19) 



e aMJ/ otherwise. 



Its action is approximately 



*7T 



2 



S in ~ — + 7r 2 p^ . (2-20) 



r 

In this background, the operator Da has a zero mode ?/>( ); 

J D^ {0) (^) = 0. (2-21) 

In the explicit notation with a for the Lorentz index and a for the SU(2) index, it behaves 
as 

at the central region and 

^}{x) ~ {2<K P )Gv af5 {x)ep a (2-23) 
at the asymptotic region, where = (<r,i), 

G FaP {x) = d^a^) aP G m= ,{x) (2-24) 

is the fermionic massless propagator, and is the two dimensional anti-symmetric matrix. 
For the path-integral of the fermionic degrees of freedom, we expand the fermion fields, 
collectively denoting them by /, as 



f(x)=J2^ n) (xH, (2-25) 

n 

f(x)=Y,V {n) \x)a{, (2-26) 
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with Grassmann coefficients a n and a n and eigenf unctions ijj^ and (p( n >; they obey eigenvalue 
equations 

-D,a,D v a v ^ = \ 2 n ^ n \ (2-27) 
-D^D^ip^ = X 2 n ^ n \ (2-28) 

and are normalized as follows; 

' d A x^\x)^ m \x) = 5 mn , f dW n) \x)<P {m) {x) = 8 mn . (2-29) 



Note is the one with A = 0. The fermionic action Sf becomes diagonalized, 

S { = £ X n d{a{, (2-30) 

n^O, f 

because the nonzero modes of (p( n > and (p^ are related one-to-one, 

¥><»> = -IfoigViW (2-31) 
The zero mode t^ -* does not have the partner and its coefficients Oq do not appear in 



Eq. fl2-30| ). The fields qi, ... I put in the integrand of Eq. (2-8) provide a proper number of 



Oq for the non-vanishing Green function against the path-integral measure; 

[dfdf\ oc daftdaff da q 3 da l . (2-32) 

We now have all the necessary ingredients to evaluate the Green function (2-8). The 
fermion fields in the integrand are now replaced with the zero mode. For the bosonic fields, 
we can just put the background instanton, Eqs. (2-14-2-17), to the leading order in g: the 
nonzero mode fluctuation around the background will give smaller contribution as long as 
n w or 77-h is not so larg e0 and zero mode fluctuations, corresponding to the translation and 
the SU(2) isospin rotation, can be taken care of collective coordinate integral as usually done 
(see the Eq. ( 2-33|) below.) The path- integral is carried out by the Gaussian approximation 



and results in the determinants from the bosonic and fermionic nonzero mode fluctuations. 
We then perform the Fourier transformation on the resulting Green function G to use the 
LSZ procedure to extract the proper vertex. The external momenta have the magnitude of 
the particles masses and are smaller than the dominant instanton size p ~ \jv. Noticing 



*' The nonzero mode fluctuations of W and 4> are of the order of g compared with the background. 
The next leading correction is obtained by replacing any two of the background with a propagator in the 
background; there are about n^y ways to pick up the two of W and ways for <j>. It is thus of the order of 
n w9 2 or n Hff 2 compared with the leading term. 
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that both the fermionic zero mode and the bosonic configuration of the instanton have the 
form of the propagator at the asymptotic region, we obtain 



r 



for the vertex, where 



dRp 2 D(p) 







n w 


n Xa f a f 










J =1 



na _ 47rV u 
Xaa = 2irpU ab € ab , 



(-2n 2 p 2 v) ni 



[2-33) 



[2-34) 
[2-35) 



and pj is the momentum of the j'-th external W boson. The integral measure dR represents 
the functional integral over the zero mode of the isospin rotation; the corresponding variables 
R ab and U ab are related by Ur b W = r a R ab . The integral over the translational zero mode has 
been done already when we Fourier-transformed the Green function to arrive at Eq. (g3p , 
and it provides the momentum conservation at the vertex. The factor p 2 ~ 5 D(p) is the 
product of the Jacobians for the bosonic collective coordinates and the determinants of the 
bosonic and fermionic nonzero modes. 

The main outcome of the above Ringwald and Espinosa analysis is that the so-obtained 
vertex (|2-33 ) is point-like independently of the momenta of external particles. Thus the 
cross section increases monotonically as ~ 

_£(2+2n w +2n H ) as the 

energy E of the process 

increases. It would eventually violate the unitarity bound if one naively extrapolate the 
result to high energy. One may think that we could predict nothing about the high energy 
behavior of the amplitude from the calculation we have done above: We have just used the 
asymptotic form of the instanton, while it is the form at the central region that determines 
the behavior at E > 1/p ~ v. This very plausible argument, however, is false. As we have 
seen, the Green function to the leading order in g is basically the product of the instanton 
configurations. The Lorentz invariance then forces its Fourier transform to depend on the 
simple momentum squares p 2 of the external particle, not on p^ ■ pj for example. Then we 
can use the asymptotic forms since we will put the on-shell condition p 2 = m 2 after all when 
we apply the LSZ procedure. 

The single instanton never disappears from the functional space no matter how large 
we take the energy in the calculation. Therefore, there must be other configurations that 
cancel it, becoming also the same order of magnitude at the energy of the unitarity violation. 
Aoyama and Kikuchi have shown that the multi-instanton configurations take this task. B 
The interaction between the instantons plays the important role in this cancelation. With 
the interactions the instantons construct "connected diagrams" for the scattering amplitude 
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and in fact they have the property that is required to recover the unitarity. Let us now turn 
to this part. 

We investigate the contribution of the n-instanton-n-anti-instanton configuration, 

n n+n 

W(x) = J2Win(x-w h Ri)+ J2 Wtnix-w^Rj), (2-36) 

i=l j=n+l 
n+n 

4>{ X ) = E ~~ (2-37) 

fc=l 

to the Green function. In Eq. (2-36), W- m represents the anti-instanton, whose form is given 
by replacing i] afll/ in Eqs. (2-14) and (2-16) with fj aiiv (fj aill> is minus of rj afJil , for \i = 4 or v = 4, 
otherwise the same.) The parameters Wi and Ri are the location and isospin- orientation of 
the i-th (anti-)instanton. We intend to treat the parameters Wi, Ri, and pi as the collective 
coordinates for the path-integral. In the strict sense, we should extract out the canonical 
variables for the fluctuations at each values of the collective coordinates and carry out first 
the functional integral over the fluctuations, and then integrate the resulting effective action 
over them. This is, however, impossible in practice. We adopt two approximations. First, 
we approximate the functional integral over the fluctuations to the simple product of the 
determinants of each (anti-)instanton. Although the configuration is not stationary in any 
sense, we do not expect this will cause a large error since the fluctuations are suppressed by g 
(the coupling constant) compared with the collective coordinates. Second, we approximate 
the action to the sum of the instanton action S- m of each (anti-)instanton and two-body 
interaction of each pair evaluated at large separation. These approximations are quite ad- 
equate if the relative distances \wi — Wj\ are larger than or m^ 1 , but if some of them 
become smaller, we will miss some part of the exact contribution. It will turn out that the 
contribution taken into account by the above procedure is sufficient to show the cancelation 
of the single instanton contribution at the energy of unitarity violation. 

The bosonic part of the interaction between the i-th instanton and j-th anti-instanton is 
defined by 

■Wi,j) = S[W(i) + W{j), 4>(i) + <f>{j)] - 2S in (2-38) 

and readily evaluated by using the fact that the each background configuration satisfies the 
equation of motion in the asymptotic region (see Ref. 0) for detail); 

d - d 

Sb:tat(»,j) = -Pv\(Pi, ^l^^ift' R j)-Q-^ G mw( W i ~ w j) 

-(2n 2 p 2 i v)(27i 2 p 2 J v)G mu (w i - Wj ), (2-39) 

where (3 is defined by (2-34) with fj substituting 77. The interaction term for a pair of 
instantons or anti-instantons is written in a similar equation in which j3 or (3 are replaced 
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correspondingly. Regarding the fermionic part, the zero mode located at each instanton 
has now overlap with the other, which also generates the interaction. Starting from the 
expansion (we suppress the suffix for the zero- modes hereafter), 

f(x) = y~] a^ l ip{x — Wj) + (nonzero modes) (2-40) 

i 

f(x) = ^ ( x ~ w j) + ( nonzero modes), (2-41) 

j 

we obtain 

= a fj a fi J dx^{j){-d^)a^{i) 

= -a h a Sl Xa a GYap{Wj - Wi)X(3a, (242) 

where x aa = 2f:pUl a e\ a comes from the fermionic zero mode ip at the anti-instanton (its form 
is Eq. ( |2-22| ) with a for a.) We have also used the fact that the zero modes satisfy Eq. ( |2-21| ) 
to go from the first to the second line. Note that the interactions (2-39) and (2-42) have the 
one-particle exchange type with the vertex exactly the same as the single instanton one for 
the relevant particle production. 

These interacting multi-instantons contributes to the B and L violating Green function. 
As an example let us see the contribution of two instantons, labeled 1 and 3, and one anti- 
instanton, 2, to the process q\qi — ► q^l (See Fig. [I| (b).) 




(a) (b) 

Fig. 1. The diagrams in the instanton expansion of the process q\q~2 — > q^L 

The corresponding Green function is 
G®( Xl ,...,x 4 ) 

= 4 / [n#]W[#ki(^i)-4^) exp {-35 in - £ S h:int (i, j) - ]T S i:ixit (hj)} ■ 
f 

(2-43) 
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To evaluate this, we expand the factor e 5int in powers of Si n t and pick up nonzero terms. 
The simplest one in G^ is 



G® 



IJ / dpid *\ dWt D( Pi )e- s ^ 



i=l 



Pi 



/ da fl dd f2 da f3 



L f 



xa qil ip(x 1 — w 1 )a g2l ip(x 2 — Wi)a qs3 i/)(x 3 — w 3 )a m ip(x A — w 3 ) 
xa 93 V 3l x(2)G F («; 2 - Wx)x{l)a l2 a ll x(2)G F (w 2 - u>i)x(l) 
xa" l2 a' ?l3 x(2)G F ( W2 - W3 M3)a 92 V 23 x(2)S F («; 2 - ™ 3 )x(3) 

i=l ^ Pi 

x (G F (m - twi)x(l)) (G F (x 2 - twi)x(l)) (G F (x 3 - tw 3 )x(3)) (G F (x 4 - w 3 )x(3)) 
x (x(2)G F ( W2 - ^i)x(l)) (x(2)G F ( W2 - ^i)x(l)) 

x ( X (2)G F (w 2 - w 3 )x(3)) (x(2)G?fK - ™ 3 )x(3)) . (244) 

This function has the same structure as the one obtained from the diagram (b) in Fig. [I] with 
the "Feynman rule", the vertices are understood as r in ( 2-33 ) or F defined for the anti- 
instanton with the proper replacements of (3 and \i an d the lines are the free propagators. 
The similar analysis can be extended for the other terms that contain any powers of S- m t- 
This way we will get the instanton expansion for the Green function. The Green function 
of a given process gets the contribution from all possible multi-instanton sectors in which 
n — n is fixed by the degree of B and L violation. Each term in the given sector is obtained 
by expanding e~ Sint in Si nt . 

We can show the terms generated by the instanton expansion is exactly the same as those 
generated by the following effective Lagrangian; B 



/ A 

+ 7(W? - d u W«Y + lm w 2 Wf + V + V 



(2-45) 



where 



V 



V 



dpdRD(p) 



5-2 



P 

dpdRD(p) 



P 



5-2 



n/(*)x 

, / 

f 



exp[-3 b 
exp[-S- u 



2^p'v)cj>(x)+P«d IM W 1 



27r 2 pM0(x)+/3« d„W2\. (2-46) 



The single instanton amplitude now turns out to be the leading term of this expansion. If one 
moves back to the Minkowski metric, this Lagrangian becomes hermitian. This indicates the 
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non-leading terms in this expansion give the contribution that cancel exactly the unitarity 
violating high energy behavior of the single instanton amplitude. We can use for example 
the K-matrix method to show this cancelation.^ [The K- matrix method provides us with 
a systematic way to pick up the diagrams in different orders of the instanton expansion so 
that their sum holds the unitarity. In the leading order of the K-matrix method, the internal 
momenta of loop diagrams are set on-shell and we can see the cancelation without suffering 
from the ultra-violet divergence that the instanton expansion has in the non-leading terms.] 
In this respect, the unitarity problem that exists in the single instanton sector is resolved. 

When returning back to the initial problem, the reliable estimate of the B and L violating 
amplitude, we note that the instanton expansion constructed so far is not sufficient. Once 
one wants to evaluate the amplitude to the non-leading orders, he encounters pathological 
ultra-violet divergences coming from the loop integrals of internal particles that connect the 
instanton vertices. The cure must come from the more rigid evaluation of the action of the 
multi-instanton configuration. Remember the approximation we adopt for the calculation of 
the interaction is adequate only for well-separated instantons, not for those with small sepa- 
rations. The interaction in the latter configuration must be important for the loop integrals 
of large momenta. We believe the proper valley method that is introduced and explained 
in the following sections is the most promising to this end. To choose the proper collective 
coordinates for the functional integral is indispensable in the multi-instanton calculation of 
the path-integral. The proper valley method has the potential to tell us the way even for 
those configurations in which a small-separated instanton-anti-instanton pair is merging into 
the vacuum. 



§3. Proper valley method 

In order to describe the proper valley method in the most general context, we will dis- 
cretize all the bosonic degrees of freedom and denote then by fa. The subscript i stands for 
coordinate labels in quantum mechanics, particle species, vector and tensor indices and any 
internal quantum numbers, as well as the space-time coordinates x M . Translation to the con- 
tinuum indices is trivial; sums over indices should be understood as integrals with suitable 
measure, the derivatives as functional derivatives.0 These variables fa are the "coordinates" 
of the functional space, which allows some intuitive understanding of the equations. The 
bosonic action is written as a function of these variables; S = S(fa). In this notation, the 
equation of motion is written as d{S = 0, where dj = d/d<fij. We assume that the metric of 

*) Inclusion of fermionic degree of freedom is done in the background of the valley configurations, as will 
be seen later in Section 4. 
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the functional space is trivial in these variables. Otherwise, the metric should be inserted in 
the following equations in a straightforward manner. 
The proper valley equationH is the following; 

DijdjS = XdiS, (3-1) 

where summation over the repeated indices are assumed and 

Dij = dtdjS. (3-2) 

Since (3-1) has a parameter A, it defines a one-dimensional trajectory in the 0-space. (The 
cases when we have a finite region with A = can be treated without any problems. This 
will be explained later in this section. For the moment, we assume A ^ 0.) The solution of 
the equation of motion apparently satisfies the proper valley equation (3-1). In this sense 
the proper valley equation is an extension of the (field) equation of motion. 

According to the proper valley equation (3-1), the parameter A in the right-hand side of 
(3-1) is one of the eigenvalues of the matrix D^. Therefore the proper valley equation requires 
that the "gradient vector" djS be parallel to the eigenvector of with the eigenvalue A. 
A question arises which eigenvalue of we should choose for (3-1). As we will show later, 
the proper valley method converts the eigenvalue A to a collective coordinate, by completely 
removing A from the Gaussian integration and introducing the valley trajectory parameter 
instead. Thus the question is which eigenvalue ought to be converted to a collective coordi- 
nate for the given theory. This really depends on the theory in question. Therefore we will 
address this issue later as we investigate various models. Here we will simply mention that 
the general guideline is to choose A to be the eigenvalue with the smallest absolute value, 
i.e., the pseudo-zero mode, or the negative eigenvalue, for which the Gaussian integration 
converges badly or diverges. In some models, the eigenvalue changes from a negative value to 
(still negative and) pseudo-zero value along the valley. In such a case, proper valley equation 
will be shown to be able to deal with both types of eigenvalues. 

The proper valley equation (3-1) can be interpreted within a framework of the variational 
method: Let us rewrite it as the following; 

9i Q (djSf - A5) = 0. (3-3) 

*) This method in the manner described here was proposed by two of the current authors, H. Aoyama and 
K. Kikuchi in Ref. [ll]). Later they found that the same equation had been written down by P. G. Silvestrov 
in Ref. |l2|), and well before that (in different contexts) in Ref. |r|). (Some very early literatures of general 
discussion of the contour-lines in geological context are found in Ref. fTi}).) Furthermore, this method is 
apparently useful in the theoretical chemistry. LLj' The authors thank Dr. Silvestrov, Dr. Quapp and S. Takagi 
at Tohoku University for letting us know of these references. 
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This allows an interpretation that the norm of the gradient vector is extremized under the 
constraint S = const, where A plays the role of the Lagrange multiplier. Since such a point 
is found at each hypersurface of constant action, the solutions of the proper valley equation 
form a line in the functional space. This is an alternative explanation for the existence of 
the valley line. In addition, we require that the norm be minimized. Since the valley with 
maximum norm is disjoint from the one with minimum norm, this can be easily done by 
choosing the appropriate valley line. We are therefore defining the valley to be the trajectory 
that is tangent to the most gentle direction. This is a plausible definition, suitable for the 
word "valley". 

What is important, however, is not its intuitive interpretation, but the evaluation of the 
relevant functional integral. It is carried out along the valley line in the following manner: 
Let us parametrize the valley line by a parameter a and denote the solutions of (3-1) as 0(a). 
The integration over a is to be carried out exactly, while for other directions the one-loop 
(or higher order) approximation is applied. We are to change the integration variables (fa) 
to a and a subspace of (fa), which is determined uniquely by the following argument: We 
first expand the action in terms of the fluctuation = — 0(a), 



The first order derivative term in the above is not equal to zero, for 0(a) is not a solution 
of the equation of motion. In such a case, the ^-expansion is no longer the loop expansion, 
as the tree contribution floods the expansion. In order to avoid this, we force this term to 
vanish by the choice of the subspace. It is most conveniently done by the Faddeev-Popov 
technique. We introduce the FP determinant A(<f>(a)) by the following identity; 



where i?j(0(a)) = diS/ y (dS) 2 is the normalized gradient vector. From the above definitions, 
we find the FP determinant to be the following; 



In the one-loop approximation, the second term in the curly bracket can be neglected and 
the remaining term is simply the cosine of the angle between the gradient vector and the 
vector tangent to the trajectory. This is simply a Jacobian factor, because the trajectory is 
not necessarily orthogonal to the chosen subspace. This situation is illustrated in Fig. 2. 



5(0) = S(<f>(a)) + diS((j>(a))fa + -D^a^fafa + ■■■. 



(3-4) 




(3-5) 





(3-6) 
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Fig. 2. A two-dimensional model of the functional space. The thin solid lines denote 
the contours of the action S. The thick solid line is the valley line. The gradient 
vector Ri and the tangent vector dcpi/da are denoted by the broken and solid 
arrows, respectively. The vertical dotted line is the subspace chosen for the 
one-loop integral. 

Let us evaluate the the path-integral for the vacuum-to-vacuum transition amplitude, 

Z = AfJl[^=e~ s , (3-7) 

using the valley trajectory, where M is a suitable normalization factor. The evaluation for 
observables is essentially the same.Ell Substituting (3-5) into (3-7), we obtain 



Z = M /da/n-^Mta) 

J j V^7T 



-j—Ri 
da 



e- sw , (3-8) 



where we have already used the one-loop expression for the FD determinant (the Jacobian) 
A(<f)(a)). We note that the above expression is apparently invariant under the general 
coordinate transformation of a. The second (functional) integration, JY[jd(f)j8((f>iRi), is 
the integration over the subspace denoted by the dotted line in Fig. 2. At the one-loop 
order, the (f>i integration yields det', the determinant of Dij restricted to the subspace. This 
subspace is orthogonal to the gradient vector, which is the eigenvector with the eigenvalue 
A according to the proper valley equation (3T). Therefore this subspace is the whole space 
less the direction of the eigenvalue A of D^. Therefore the resulting determinant is simply 
the ordinary determinant less the eigenvalue A; 

det'D = ^ (3-9) 
A 



*) The only possible exceptions are the cases then we have "large observables" ; if for example we have 
a large number of fields, which contributions could override the dominations by the exponentiated action. 
This requires further development of the theory at hand. 
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This can be shown by evaluating (3-8) by exponentiating the delta-function and integrating 
over the fluctuations 6 as follows; 



o v 

00 dk 



z=n da /n^= 

= AT J da J 
= M J da 



dk 

2n J-oo 2tt 
1 



,ik<piRi 



da 



Ri 



-s(<t>) 



, exp i --RD^Rk 2 

2tt y/MD V 2 



da 



R, 



2n(RD- l R) detD 



(3-10) 



Since the proper valley equation (3-1) yields D 1 R = R/X and the vector R is normalized, 
we obtain 

(RD- l R) det D = - detD{= det'D), (3-11) 
A 

as in (3-9). We thus obtain the following expression at the one-loop order; 



Z = N J da 



-j—Ri 

da 



V2n det' D 



S(<Ka)) 



(3-12) 



In this sense, the proper valley method converts the eigenvalue to the collective coordinate, 
exactly. As we briefly mentioned at the beginning of this section, this exact conversion is 
quite ideal for the actual calculations: In physical situation one often suffers from a negative, 
zero, or positive but very small eigenvalue, which renders the ordinary Gaussian integration 
meaningless or unreliable. The proper valley method saves this situation by converting the 
unwanted eigenvalue to the collective coordinate, for the factor det' D is exactly free from 
this eigenvalue. 

When there are solutions of the equation of motion with a zero mode, as in the models 
with ordinary instantons, the proper valley method reduces to the usual collective coordinate 
method. This is seen as follows. Let us denote the solution of the equation of motion by 

0soi; 

diS{<f> sol ) = 0. (3-13) 

This satisfies the proper valley equation (3-1) at the same time. The normalized gradient 
vector Ri is apparently singular due to ( |3-13| ). However, Ri is well-defined by a limiting 
procedure A — > 0. In order to show this, we denote the original action by and we add a 
term XS^ to the original action, so that so i is no longer a solution of equation of motion 
for the theory defined by the total action S = + XS^ 1 '. (The parameter A is an arbitrary 
small quantity and can be any function of <p that satisfy the above requirement.) The 
solution of the proper valley equation for 5* can be expanded in terms of A around (j> so \. 
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For this purpose, We introduce a set of auxiliary variables Fi and rewrite the proper valley 
equation (3-1) as follows; 



DijFj = XFi, (3-14) 
d i S = F i . (3-15) 



The equation ( 3-13 ) guarantees that at A = 0, Fi = 0. Thus we use the following expansion; 



4> = so i + A0 {1) + • • ■ , 
F = AF (1) + • • • . 

Substituting the above expansion to (3-14) and (3-15), we find that the zeroth order equation 
in A is simply fl3T3|) , and that the first order equations £1X6 £IS the follows; 



D^fV = 0, (3-16) 
^% sol )+£gVf =if\ (3-17) 



where = didjS^((f) so i). From the equation (3-16), we find that F^ is proportional to 
the zero-mode of the original solution so i. The proportionality constant is determined by 
multiplying F^ on (3- 17); 

F/ 1) ^S«(0 sol )=F/ 1) F/ 1) , (3-18) 

where we have used the fact that Dij is symmetric and (3-16). The quantity (jyp is determined 
by (3-17) under the above condition. Since the normalized gradient vector Ri in question is 
the normalized F vector, i.e., Ri = Fi/\F\, it is, in the limit A = 0, simply the the normalized 
zero-mode of the solution of the equation of motion. Therefore, we conclude that the FP 
procedure (3-6) is equivalent to the usual collective coordinate method, as promised. 

The above analysis shows that the proper valley method can be considered as an gener- 
alization of the usual collective coordinate method. An added bonus of this property is that 
the resulting det'-D is quite easy to calculate; we simply calculate the whole determinant 
and divide it by the smallest eigenvalue. If several unwanted eigenvalues exist, the proper 
valley method can be extended straightforwardly. The equation should then specify that the 
gradient vector diS lie in the subspace spanned by the unwanted eigenvalues. This leads to 
a multi-dimensional valley with the advantages noted above. 

There is another valley method, called "streamline method" 0), which has been exten- 
sively used in the literatures. It proposes to trace the steepest descent line starting from a 
region of larger action. Namely, it uses the equation 

^ = f^S, (3-19) 
da 
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where f{a) is an arbitrary function of a (there is a general reparametrization invariance of 
a in this formalism). This valley line coordinate can be introduced in the same manner as 
before using the FP technique. This makes the Jacobian trivial at the one-loop order. This 
streamline method, however, suffers from many difficulties unknown to the proper valley 
method: The most serious one is the fact that its subspace for the Gaussian integration 
has no relation to the eigenvalues of Dij whatsoever. The relation ( BvTT ) does not hold; the 



determinant is not guaranteed to be free from the unwanted eigenvalue. Another problem 
is that it is a flow equation in the functional space: Once given a point in the functional 
space, it defines other valley configurations in the infinitesimally small neighborhood of that 
point. In other words, it is not a local definition in the functional space. As such, it does not 
define a solution by itself. One has to define the starting point to start the valley. Therefore 
the construction of the configurations on the valley trajectory is quite difficult. When the 
starting configuration is known (like a pair of instanton and anti-instanton separated by 
infinite distance), this might be said to be a technical difficulty In this method, however, 
the valley can be traced only from the above; it suffers from instability if the valley is traced 
from the bottom of the valley. An elementary analysis shows that the streamline most often 
climbs up the side-walls of the valley, not the bottom line. This problem becomes fatal when 
one is dealing with the theories with zero-size instantons, like the standard model noted 
in the previous section, in which the higher starting point is not known at all. For more 
detailed comparison of these methods and other features of the proper valley method, we 
refer the readers to Ref. |TT| ). 

There is yet one another formalism that deals with this kind of situation. It is the 
"constrained (instanton) method", 0* explained briefly in the context of the standard model 
in Section 2. This method is actually a valley method, in the sense that we obtain a line of 
configuration parametrized by the constraint parameter. In the case mentioned above, this 
method allows the definition of configurations, are very close to the zero-size instanton. It 
has a problem that there is no prescription for the choice of the constraint that guarantees 
the effective evaluation of the path-integral. In general one knows some way to choose the 
constraint so that it allows desirable configurations. But this is only a prescription for 
the constraint not to be meaningless. Among all the meaningful constraints which allows 
effective evaluation is the question unanswered. 



Actually there is a subtle relation between the proper valley method and the constraint 
method. Let us elaborate upon this point for clarification. The equation (3-3) can be 
interpreted as a constrained equation in a different manner: By dividing with A, it can be 
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viewed as extremizing the action S under the constraint functional; 



S, 



const. 



(3-20) 



Therefore, the proper valley method can be interpreted as a constraint method, but with 
the specific choice of the constraint G3-20|) . However, with all the geometrical properties 
and the special features, especially in relation with the Gaussian integrations, which are all 
unique to the proper valley method, we consider it is best to distinguish it from the generic 
constraint method without any special consideration to the constraint. Hereafter in this 
paper we simply call the latter the constraint method. 

While the detailed comparison of the configurations obtained by the constraint method 
and the ones by the proper valley method in specific field theoretical models are done in 
Section 4, here we discuss a toy two-dimensional functional space to demonstrate the power 
of the proper valley method. The two degrees of freedom in the model is denoted as 0j 
(i = 1,2), and the action is given by the following, 



5(0i,0 2 ) = 4((^i + ^) 2 + 5(0?-0 2 
9 v 



(3-2i; 



This is constructed by distorting a simple parabolic potential so that the valley trajectory 
is not trivial. Its contours (S = constant lines) are plotted as thin lines in Fig. 3. 

02 



0.5 




0.5 0i 

Fig. 3. The valley instanton and the constrained instanton in the toy model. The 
solid line denotes the valley instanton and the dashed line shows the constrained 
instanton. The dots show the positions of the saddle points of the integrand of 
0(p,q). 

The proper valley equation is now a simple algebraic equation, which can be solved 
numerically. Alternatively, we could reparametrize the functional space by (r, 9) defined by 
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the following; 

2 + 2 = r C os 9, <f>\ - 02 = -^rsin#. (3-22) 

Under this parametrization S = r 2 /g 2 . Therefore following the variational interpretation of 
the proper valley equation, we can minimize the square of the norm of the gradient vector, 
(diS) 2 , as a function of 9. In Fig. 3, the thick solid line denotes the proper valley trajectory. 

As an analogue of the constrained instanton formalism in this toy model, we introduce 
a simple constraint. Near the origin the valley trajectory extends to the 0i direction. Any 
reasonable constraint has to reproduce this property. Therefore, as a simple example for the 
constraint, we choose the following; 

0i = const. (3-23) 

The solution of the constraint method is plotted in Fig. 3 by the dashed line. 

In Fig. 3, it is apparent that the proper valley trajectory goes through the region of 
importance, while the constrained trajectory does not. This property can be displayed more 
explicitly by considering the analogue of the physical observables 0(p,q) defined by the 
following; 

0(p,q) = #! / rf0 2 0?01 e ~ s ^\ (3-24) 

JO J-oo 

Following the prescription given above, we have carried out the numerical evaluation of 
0(p,q) exactly and by the Gaussian ("one-loop") approximation around the proper valley 
trajectory and the constrained trajectory. Tabled gives the ratio of the one-loop values over 
the exact value. 

As is seen in the table, the proper valley method gives a better approximation than the 
constraint method in this range of (p, q) consistently. This could be explained in the following 
manner; one could calculate the position of the maximum of the integrand in (3-24), or the 
following "effective action"; 

£(™)(0) = 5(0) -p log 0! - q log 2 . (3-25) 

The peak positions calculated in this manner are denoted in Fig. 3 by dots. The point 
a in the left bottom is for (p,q) = (2,2), b for (2,12), c for (12,2), and d for (12,12). 
The dots are distributed around the proper valley trajectory. This explains the fact that the 
proper valley trajectory gives the better approximation over the constrained trajectory. Only 
when the power of one of the observable gets high enough as in the point c, the constraint 
method starts to give a little better values. However, this is exactly the range in which the 
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ordinary instanton calculation starts to fail; For the n-point function with n > 0(1/ g 2 ), it 
is well known that one needs to take into account the effect of the external particles on the 
instanton itself. The power of the valley method lies in the fact that even in this range the 
calculation is done with good accuracy for generic type of observables, like the (12,12) case. 
Only when the observable is very special, like the (12,2) case, the constraint method tuned 
to that operator may give reasonable estimates. 
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0.827 


0.679 


0.585 


0.520 


0.474 


0.439 




0.545 


0.248 


0.110 


0.049 


0.022 


0.010 


4 


0.997 


0.870 


0.767 


0.689 


0.629 


0.581 




0.667 


0.324 


0.149 


0.068 


0.030 


0.014 


6 


1.061 


0.986 


0.900 


0.826 


0.764 


0.712 




0.737 


0.384 


0.185 


0.086 


0.040 


0.018 


8 


1.064 


1.046 


0.991 


0.931 


0.876 


0.826 




0.783 


0.432 


0.218 


0.105 


0.050 


0.023 


10 


1.029 


1.063 


1.041 


1.003 


0.960 


0.918 




0.815 


0.473 


0.248 


0.123 


0.060 


0.028 


12 


0.971 


1.046 


1.058 


1.043 


1.017 


0.985 




0.839 


0.508 


0.276 


0.141 


0.070 


0.034 



Table I. The ratio of the 0(p,q) estimated by proper valley method (the upper 
column) and the constraint method (the lower column) over the exact value for 
5 = 0.2. 

Let us note a tricky point in comparing the valley instanton and the constrained instanton. 
If we compare them with the same parameter values, the constrained instanton has smaller 
action than the valley instanton. This is not a contradiction, nor it means the constrained 
instanton has a larger contribution: This becomes clear by considering the dotted vertical line 
in Fig. 3. The point x and y are the valley and the constrained configurations, respectively, 
which have the same value of 4>i- The configuration y has a larger action than x, by the 
definition of the constraint method. However, as we have seen above, this does not mean 
that the constrained trajectory gives a better approximation. (One way to illustrate this 
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point is that if we compare the configurations at the same distance, / da\{d<f>i/ da)Ri\, from 
the origin, the valley configuration has the smaller action.) The same situation will be seen 
in the results of the following sections; when we compare the instantons at the same value 
of the constraint parameter, the constrained instanton has a smaller action than that of the 
valley instanton. This is a red herring and has nothing to do with the relevancy of the valley 
instanton. 

In this analysis we took a very simple constraint (3-23). Alternatively we could take 
some ad-hoc constraint as long as any of these yield "finite-size" instantons, i.e., points 
away from the origin. Most of these constraints yield essentially similar results. Only when 
the trajectory goes through the dotted area in the Fig. 2, one can obtain good quantitative 
results. Such a trajectory, however, is destined to be very close to the proper valley tra- 
jectory. Therefore the proper valley method, defined without any room for adjustment and 
guaranteed to give good results is superior to the constrained method. 

Finally, we note that in the actual calculation of the solutions in the continuum space- 
time the proper valley equation with auxiliary variable F, the equations (3-14) and (3-15), 
is most convenient. Let us denote all the real bosonic fields in the theory by <f> a (x), where 
we denote all the space-time coordinates by x and the rest of the indices by a. The proper 
valley equation is as follows; 



This is a set of second-order differential equation, which can be analyzed by the conven- 
tional methods. The solution of the ordinary field equation of motion is a solution of the 
above equation with F a (x) =0. In other words, F a (x) specifies where and how much the 
valley configuration deviates from the solution of equation of motion. This property is useful 
for the qualitative discussion of the properties of the valley configurations. The analysis in 
the following sections will be carried out in this auxiliary field formalism. 



4.1. Valley instanton in the asymmetric double-well potential 

In this subsection we apply the proper valley method to an one-dimensional quantum 
mechanical model. Our aim is to resolve a paradoxic that arises if we naively apply the 
semi-classical approximation to the potentials, V(<p), depicted in Fig. 4 and 5. Let us start 




(3-26) 



F a (x) 



5S 



8<j) a {x) 



4. Application of the proper valley method 
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with reviewing the path-integral formalism in the quantum mechanics and then clarify what 
is the paradox we are concerned with. 




For a given Hamiltonian with the dynamical variable <fi and its conjugate momentum p, 

1 



H= -P 2 + Vtf), 



(4-1) 



the transition amplitude from the initial value <fi = qy at r = to the final one q 2 at T in 
imaginary time formalism is evaluated by the path-integral 



(q2\e- HT/h \ qi ) = M Jv<t>e 



sm/n 



where 



S[<j>] 



dr 



In Eq. (|4-2|) , <j){r) is parameterized as 

n 

where 0b is a background that satisfies 

0b(O) = qi, B (T) = q 2 , 



(4-2) 



(4-3) 



(4-4) 



(4-5) 



and are elements of a normalized complete set in the functional space for the interval 
[0, T}; they satisfy 



and 



(p<ri(0) = <pW(T) = 0. 



(4-6) 
(4-7) 
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The integral measure is then defined in terms of the coefficients a n by 

(4- 



da 1 x 2 1/2 

n V2liTl 



1 fi(T' 



2i,hT *i, V T 

The infinite products over m and n are understood to be combined properly in actual calcu- 
lations so that the final expression becomes a convergent product. The information on the 
model will be given by calculating Eq. ( J4-2|) for the various values of q±, q 2 , and T. We are 
interested in the energy eigenvalues E of the lowest-lying states; we calculate Eq. ( |4-2|) with 
the boundary condition 0(0) = 0(T) = and with sufficiently large T (compared with the 
inverse of the excitation energies of the other states.) The energy eigenvalues are extracted 
from the T dependence of the transition amplitude, (0\e~ HT ^ h \0) ~ e~ ET l n . 

The semi-classical approximation for the path-integral Eq. (|4-2|) is carried out by the 
following: use the classical solutions c i 

SS 



as 4>b and expand the action around it, 



(4-9) 



S[<j>] = S[<j> d ] + ~ J J drdr'D(r, r')^(r)^(r'), (4-10) 

where D(t,t') = 5 2 S/S4>(t)S(P(t i ) and ip = — 4>c\] use the eigenfunctions of D for the 
complete set {v?^} and approximate the functional integral by Gaussian, 

(0|e^ T |0)^E^ l] 7^= (441) 

/detD[0 d ] 



where "det" denote the product of all the eigenvalues of D (we will set % = 1 hereafter.) The 
zero modes will be treated by the collective coordinate method (see for example Ref. [|)) if 
D has some. 

For the potentials in Fig. 4 and 5 we have the bounce solution, that starts at <fi = and 
bounce at 0esc back to (ft = for the large time interval, T ^> u^ 1 (u + is the frequency 
at = 0.) There is a subtlety in the semi-classical approximation at the bounce: it has 
one negative eigenvalue mode and the Gaussian integral along this direction becomes ill- 
defined. The way we usually adopt to circumvent this is to assume an appropriate analytic 
continuation of the contour for the ill-defined Gaussian integral; The bounce then gives an 
imaginary part for the energy eig envalues!)'@)>@> for both potentials in Fig. 4 and 5. 

Although the exact form of the analytic continuation is not so clear, the complex energy 
for the potential in Fig. 4 has a sound interpretation: if one restricts the wave functions 
to have only the outgoing component at ^> 0esc ; the hermiticity of the Hamiltonian 
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is violated and the energy eigenvalues become complex. This imaginary part reflects the 
instability of the localized wave packet at ~ 0. 

For the potential in Fig. 5, however, the imaginary part is obviously a wrong answer. 
We can only take wave functions with decaying exponential at 0o an d the hermiticity 
of the Hamiltonian cannot be violated. There is no room for the energy eigenvalue to be 
complex. This is the paradox we intend to solve. 

We use the proper valley method to see how the action behaves for the path-integral 
related to the negative mode at the bounce and to make the corresponding path-integral 
well-defined. We construct "valley instanton", which should replace the bounce solution. 
Interestingly, it has a zero mode and this expedites the calculation of its determinant and 
Jacobian as was the case of the instanton. It is a well-localized configuration with respect to 
the imaginary time and their dilute-gas sum generates the reasonable energy shift instead of 
the decay rate. We will also show that it converges analytically to the instanton in the limit 
of e — > and all the results reproduce the well-known instanton results. 

Specifically we parameterize the potential in Fig. 5 as 

V(<f>) = (1 - g<pf - e(4/0 3 - 3^ 4 4 ), (442) 

where the coupling constants g and e are positive. The potential ( |4-12| ) has a local minimum, 
V(0) = and a global minimum at O = l/g, where V(l/g) = — e. (Fig. 5 is plotted for 
g = 0.3 and e = 0.25.) The potential V(4>) in Eq. Q4-12Q is a canonical form of quartic 
potentials with two minima, since any such potential can be cast into this form by suitable 
linear transformations on and r. In this sense the following analysis is quite a general one. 
(Especially, the above form is related to that of Ref. |20|) by a simple reparametrization.) In 
the following, we consider the cases with small coupling g <C 1, but not necessarily small e. 
The proper valley configurations are given by the equations 

-d 2 r( p + V'(4>) =F, (4-13) 
(-fi£ + V"{<t>)) F = XF. (4-14) 

as has been explained in the section 3. The bounce satisfies this equation trivially since 
it has F = 0. We take A in the infinitesimal neighborhood of the bounce as its negative 
eigenvalue. Then the valley trajectory extends to the direction of the negative mode from 
the bounce. By changing the parameter A, we will obtain a trajectory in the functional 
space. The general solutions can be parameterized by the eigenvalue A, or any arbitrary 
function a of A. We denote the valley trajectory by 0(a). We consider the contribution of 
such a configuration to the transition amplitude (0|e~ HT |0). By the result of Eq. (3-12), this 



28 



amplitude is given as 



(0\e- HT \0)=Af Ida I J a < 
J V27rdet D 



(4-15) 



to the leading order in g, where det' D is the usual determinant less the eigenvalue A in the 
proper valley equation (4-14) and 



J a — 



_/ 



d(f)(a) dS[(f)(a)] 
aT — : F 



da 



da 



drF 2 



(4-16) 



drF 2 



The factor J a is the Jacobian for this change of the integration variable which is explained 
at the previous section. 

We have carried out the numerical analysis and obtained the solutions of the proper 
valley equation plotted in Fig. 6. The solid line, a, is the bounce solution of the equation of 
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Fig. 6. Valley bounce solutions (4>(t), F(t)) of the proper valley equations. Center of all of the 
configurations are chosen to be at the origin, r = 0, around which the solutions are symmetric. 
The solid line a in the upper figure is the usual bounce solution, which has F(t) = 0. The 
other lines, b-e, are unique to the proper valley equations. 



motion. The rest do not satisfy the equation of motion. We call these solutions (including 
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the bounce solution) "valley bounce". The values of the action, S, and the eigenvalue, A, for 
the valley bounces are plotted in Fig. 7. The bounce solution lies at the top of the line of the 



s 




Fig. 7. The values of the action, S, and the eigenvalue, A, of the valley bounces. The peak of the 
action is given by the bounce solution, the line a in Fig. 6. The points corresponding to the 
valley bounces in Fig. 6 are plotted with circle. The solid lines are drawn as the guide for eyes. 
The valley parameter is chosen to be |0| = / drcf). 

action in Fig. 7, corresponding to the fact that it has a negative eigenvalue. This negative 
eigenvalue can be read from the corresponding point of the plot of A in the lower half of 
Fig. 7. The most notable feature is that the large size valley bounces have a clean interior, 
where = O an d F — 0. (This property is shared by the higher dimensional configurations, 
i.e., the valley bubbles.) The effect of this is apparent in the behavior of the action in Fig. 7; 
the action decreases linearly with large \<f>\, which is proportional to the size of the valley 
bounce. 

In the large valley bounces, c, d, and e in Fig. 6, we notice that the shape of the wall, 
i.e., the transition region from ~ O to ~ 0, is almost identical to each other. That is, 
they overlap with each other very well when translated in r. Thus these large valley bounces 
can be approximated by simply connecting the walls by a flat region, = <f> , at various 
separation. The shape of the wall can be most readily identified when the size (|0| = / drcfi) 
of the valley bounce becomes oo. In this limit, the wall is simply a localized transition 
from = to = 0o (or vise- versa). Such a solution is an analogue of the instanton (or 
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anti-instanton). The difference is that now it is not a solution of equation of motion, but 
is a solution of the proper valley equation. This kind of configuration is called the "valley 



The notable information from this numerical analysis is that the path-integral for the 
negative mode direction at the bounce can be translated eventually into the collective co- 
ordinate integral with respect to the relative distance between the valley instanton and the 
valley anti-instanton. This integral indeed gets the large contribution from the large sepa- 
ration and is divergent for T — > oo, but we can handle it anyway; we will be even able to 
carry out the dilute-gas summation for the multi valley instantons. Before doing so, we look 
into the details of the valley instanton which is needed for the summation. 

Since we define the valley instantons at the large-size limit, \<fi\ — ► oo, of the valley bounce, 
the plot of the eigenvalue in Fig. 6 implies that the eigenvalue A of the valley instanton is 
exactly zero. This is not a trivial property. A solution of equation of motion is guaranteed to 
have zero modes corresponding to its symmetry transformation, such as a time translation. 
Arbitrary background configurations do not have this property in general. However, we 
can prove the existence of the zero mode as follows: Take a derivative of Eq. (4-13) with 
respect to r, multiply F, and integrate over r. After partial integrations (which surface 
terms vanish), we then find that 



The integral in the left-hand side is generally non-zero due to the boundary conditions of 
the valley instantons. This can be seen from the behavior of and F in the walls in Fig. 6. 
Since the right-hand side is zero, we find that A = 0. (For the valley bounces, the integral 
in the left-hand side is zero since 0(r) and F(t) are even functions. Therefore, A 7^ is 
allowed.) 

We have carried out the numerical investigation of the valley instanton with A = and 
have successfully obtained the solutions, which have turned out to be almost identical to the 
wall regions in Fig. 6. 

Although we know of no exact analytical expression of the valley instanton, it can be 
constructed analytically for small g 2 in a manner used in the construction of the constrained 
instanton0) as well as other types of the valley instantons. 0* Consider the valley instanton 
in t G [— T/2,T/2] (T ^> 1). We define its central coordinate to be at the origin, r = 0, by 
0(0) = l/2g. The naive perturbation in eg 2 yields the following perturbative valley instanton 
solution; 



(anti-) instanton." 




(4-17) 



= ^ + 3eA0o + O((es 2 ) 2 ), 
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(4-18) 

F = -6eg 2 ^ + 36^^ + 0((q? 2 ) 3 ), 
where 0q is the instanton solution for e = 0; 

K = --r^ — • ( 4 ' 19 ) 
gl + e~ T v 1 

From Eq.(4-18), it is apparent that this naive perturbation is valid only in the region close 

to the instanton center, |r| <C l/(eg 2 ). On the other hand, in the asymptotic regions, 

t — > ±oo, we linearize the proper valley equation and find the general solutions valid for 

\t\ ^> 1. Coefficients of the general solutions are fixed by matching it with the inner solution 

Eq.(4-18) in the intermediate region, 1 <C \r\ <C l/(eg 2 ). We have found that this procedure 

can be done consistently. The resulting asymptotic behaviors are for r — > +oo; 



■oj_r 



(4-20) 



F ~ -Qege- u -\ 
where u 2 = V"(4>q) = 1 + 12e5f 2 , and for r — > — oo; 

<f) ~ - fl + 3e^ 2 r) e T , 

(4-21) 

F ~ — 6e^e T . 

This way, the valley instanton is constructed in all regions of r for small e. 

The action of the valley instanton is given by S 1 = 1/Qg 2 + e(— T/2 + 1/2) + 0(e 2 ). This 
action is divided to the volume part and the remaining (proper) part as S 1 = —eT/2 + S 1 . 
From the construction above, we find that S 1 = l/6g 2 +e/2+- • •. However, there is a subtlety 
on this point: In the following we integrate over the position coordinate of the instantons 
and anti-instantons in the dilute-gas approximation. These coordinates are originally the 
valley parameters (as) of the valley bounces (and their central coordinates). The 0(e) term 
in S 1 depends on the definition of these valley parameters, since the definition of the volume 
(T) is affected by it. Therefore, careful study of the small valley bounces is needed to fix 
this term. This term, however, has only the non-leading contribution. Therefore, we will 
not pursue this problem any further here. 

The Jacobian for the instanton position is given by, 

J 1 = , 6 . (4-22) 
drF 2 
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Since the contribution to the integration is dominated by the central region, the leading 
term of Eq.( [4-22| ) for small e is evaluated by the use of Eq.(4-18). The result is that J = 
1 /v / % 2 (l + 0(eg 2 )). The first term is the instanton action for e = 0. Therefore, in the limit 
e — > 0, the Jacobian of the valley instanton reduces to that of the ordinary instanton. 

The determinant, det'D, can be calculated by extending the Coleman's method, & in 
spite of the fact that the valley instanton is not the solution of the equation of motion. This 
is due to the fact that the valley instanton possesses the exact zero mode F(t). We define the 
asymptotic coefficients F± by F(r) ~ F±e Tw±T , where for the sake of notation we introduced 
uo\ = V"(0) = 1. After some calculation, we find that the ratio of the determinants for the 
valley instanton located at r (G [— T/2, T/2]) is given by the following; 



det(-d 2 + oo 2 + ) 2uj + uj-F + F_ 



oo 



det'M? + V(f)) _ Ke ,„__„ +)W2 _ T0 , K s 1 , - pi (4 . 23) 



oc 



The exponential factor represents the contribution of the difference between the zero point 
energies at the two local minima. The factor k is the 'proper' instanton contribution. From 
Eq.(4-20) and Eq.(4-21), we find that k reduces to the ordinary instanton determinant for 
e -> 0. 

We combine all factors and calculate the amplitude using the "dilute-gas" valley instanton 
approximation. We find the transition amplitude to be the following; 

where we have normalized the amplitude to that of the free Hamiltonian H = (1/2) [p 2 + 
u/j.0 2 ]; n is the number of the valley instanton pairs, and the factor a is the product 
of the proper contributions of the action, the determinant ratio and the Jacobian; a = 
( J 1 /\ / 2nK)e~ s ' . The actual integrations over the positions of the valley instantons are in 
the factors I n ; 

!1, for n — 0, 

/ dr 2n / dT 2n -i- / dT 1 e eiT2n ~ T2n - 1+ - +T2 - T1 \ for n > 1, 
Jo Jo Jo 

where zero-energy contributions of the determinants are absorbed in e by e = e — (uo_ — uo + ) /2. 

The infinite series can be summed by the use of the generating function method:© From 
Eq.(4-25), we find that the following differential equation is satisfied by Z(T); 

Z(T)" - eZ(T)' - a 2 Z(T) = 0. (4-26) 
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Also, Z(0) = 1, and Z'(0) = 0. Therefore, we find that 



h, p k—T h, p k+T 

Z(T) = - , (4-27) 



where, 



fc ± E--±J i + a 2 (4-28) 



Thus we find that the energies of the two lowest states, E± is given by 



^f + ^-Hf + ^. (4-29) 

Furthermore, from the coefficients of the respective exponents of Eq. (|4-27 ), we find that 
|(0 = 0|i?-|-)| 2 = ±k±/(k + — k-). Note that there appear no fake imaginary parts in the 
energy spectrum. 

We also examine the validity of the dilute-gas approximation. The mean size of the 
bounce, which is the mean distance, R, between the instanton and the anti-instanton located 
at the right of the instanton, can be obtained as the expectation value, (r 2 — 7~i) in the 
amplitude, Eq. ( f4-24| ) and Eq.(4-25). (Any other (r 2m — T 2m -i) with integer m would yield 
the same result for T — > oo.) Using the same generating function method as above, we 
obtain the following; 

2 

Similarly, the mean distance, d, between the anti-instanton and the instanton located at the 
right of the anti-instanton (r 3 — r 2 ) is given by, 




(4-30) 



For e>a, d~ 1/e. On the other hand, the thickness of the instanton is 0(1/^) for e 3> 1, 
which can be seen by a simple scaling argument on the proper valley equations, Eq.(4-13) 
and Eq.(4-14). This means that the dilute-gas approximation is valid for e < 1. 

The result Q4-29|) has now a simple explanation. We observe that E± are equal to the 
eigenvalues of the following matrix; 

H= 2 V • (4-32) 

Furthermore, the weights of the state localized in the left well agree with Eq.( |4-27j ): Denoting 
the eigenvectors of Eq. ( |4-32j ) V± with eigenvalues E± , 



iWraWi-^- (4 - 33) 
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The energy spectrum we obtain is the same as the two-level system made of the perturbative 
ground state at = and <fi = (f>o, with the tunneling matrix element a. This simple picture 
is correct if the higher excited states have large energy gap from the two, and the condition 
is exactly the one obtained from the validity of the dilute-gas summation, e < 1. 

Finally let us check the result (|4-29|) for a few cases. For e <C a, the energy spectrum 
Eq. (|4-29|) gives the following; 

E ± = U++ A U - ±a-i + o(^). (4-34) 



In the limit e — > 0, this result reduces to the well-known instanton result for the degenerate 
case. In addition, we have the average of the perturbative zero-point energies of the left- well 
uj + /2 and the right-well uj-/2 — e. Since the wave function is distributed evenly at the zeroth 
order of the e expansion, this is the correct formula for the two lowest energy eigenvalues. 
Since the perturbative contribution to the energy splitting, 5E = E + — E_, is of order eg 4 , 
purist a> would retain our result, 8E = 2a. For a < e < 1, Eq.( |4-29|) leads to, 



even a 



£+ = f + T + bj' e - = -' + --- + {-J)- (435) 

The lower energy, E_ corresponds to the state almost localized in the right well; this is 
also an expected result since the tunneling transition is small compared with the energy 
difference. 

In summary, we have developed the method to obtain the correct energy shift in the 
asymmetric double-well potential by applying the proper valley method. In view of this 
development, the ordinary calculation of the imaginary part of the energy in the unstable 
cases (Fig. 4) needs to be examined under the new light. 

4.2. Valley bubble 

In this subsection, we study another example of quantum metastability problem, the false 
vacuum decay. It also can be treated in the imaginary-time path-integral formalism using 
the bounce solution. EH^ Il^ 

Due to the existence of a negative-eigenvalue fluctuation-mode around the bounce solu- 
tion, the contour of the Gaussian integration through the bounce has to be deformed to yield 
the imaginary part of the energy level. Thus the decay rate of the false vacuum is obtained. 

This bounce belongs to a valley of the action in the whole functional space. © This 
situation is analogous to that of the tunneling processes via instanton, such as the baryon 
number violation process in the standard model explained in section 2. 
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In either case, when the initial condition is such that the state is the local minimum, the 
solutions of the equations of motion (bounce or instanton) dominate the relevant imaginary- 
time path-integral. However, when the initial state is of higher energy, different configura- 
tions could dominate. For the quantum tunneling between degenerate vacua, the deformed 
instanton and anti-instanton pair is known to play that role.© These configurations gen- 
erally belong to the valley of action, since they form a line of relatively small actions. 0* 
This is obtained by separation of a collective coordinate that corresponds to small, zero or 
negative eigenvalue, which is dangerous for the Gaussian integration. On the other hand, not 
much is known for quantum decay of metastable states: In the false vacuum, decay occurs 
through a vacuum bubble.© When the energy is high, we expect that modified vacuum 
bubbles dominate the path- integral, just like instanton anti-instanton with shorter separa- 
tion dominates as energy increases. In fact, one may estimate the n-point Green function 
in a thin-wall bubble background as follows; Far from the bubble at the origin, the scalar 
field 4>(x) behaves as ~ e R ~ p for p = \x\ — > oo, where we denote the radius of the bubble 
by R. Thus the n-point function has the integrand, e nR ~ s . Minimizing this expression with 
respect to R, we find that the radius shifts by —0(en/Sf), where Si is the surface energy, 
and e is the energy density difference. This leads to a smaller bubble, whose mid-section is 
a three-dimensional bubble with energy of 0(n). This argument demonstrates that in fact 
the induced decay of the false vacuum is mediated by smaller bubbles under the presence of 
low energy incoming particles. However, this is quite a crude estimate. Accurate analysis 
should be carried out on the basis of the imaginary-time path-integral formalism. The start- 
ing point is then to clarify the structure of those modified bubbles. This series of bubbles 
forms a valley, on which the bounce is the saddle point. In this subsection, we carry out the 
analysis of this valley, employing the proper valley method. We consider a quantum field 
theory of a scalar field 4>(x) in 3+1 dimensional space-time with an action S, which we shall 
specify later. The proper valley equation (3-26) can be obtained by varying the following 
total action, S to t, 



'tot — 




(4-36) 



where 



5, 




(4-38) 



(4-37) 



We consider the following action, 




(4-39) 
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Since the above potential is same as Eg. ( |4T2| ) of g = 1, it has the false minima at <fi = and 
the true minima at <fi = 1, regardless of the value of the parameter e(> 0). The energy density 
of the true vacuum is — e, even for large e. In this sense, this action defines a convenient 
model for study of thick-wall bubbles as well as thin-wall ones. 




Fig. 8. Potential V(cf>) for e = 0.25. 

As the bounce solution is spherically symmetric, we shall confine ourselves to the study 
of spherically symmetric configurations, (j){x) = <p(p) and F(x) = F(p), where p 
The proper valley equations lead to the following; 

0" + ~<t>' " -77 + F = 0, 

p d(p 



(440) 



F" + -F' - -jj^F + XF 
p d<p z 



0. 



(4-41) 



where we denote the derivatives with respect to p by primes. Just as in Coleman's treatise 
of the bounce equation, the above equations can be thought as the set of Minkowskian 
equations of motion of a particle in two dimensional space (0, F) at "time" p. The linear 
differential terms of <fi and F act as friction terms. (Note that there are no friction term in 
Eq.(4T3) and Eq.(4-14).) The other terms are space-dependent forces. The big difference is 
that now this force is not conservative. Therefore, no simple energy argument is possible. 

The proper valley equations Eq.(4-40) and Eq.(4-41) require four boundary conditions. 
They are as follows: For the solution to be regular at the origin p = 0, we require boundary 
conditions </>'(0) = F'(0) = 0. The outside of the bubble has to be the false vacuum, so 
0(oo) = F(oo) = have to be satisfied. 

In solving Eq.(4-40) and Eq.(4-41) numerically, we have chosen to start at the origin 
with the boundary conditions <fi'(0) = F'(0) = and adjust 0(0) and F(0) so that 0(oo) = 
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.F(oo) = are (approximately) satisfied. We have done this calculation for e = 0.25. The 
potential for this value of e is given in Fig. |8|. Since the depth of the true vacuum is much 
more than the height of the potential barrier, we expect that the bounce solution is a bubble 
with a thick wall. 

In Fig. 9, we show the numerical solutions and Fig. |1TJ is the relation between smallest 
eigenvalue A of D and the starting point 0(0). We observe the following in these figures: 



4> 
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Fig. 9. Shapes of the solutions of the proper valley equation Eq.(4-40) and Eq.(4-41). The eigen- 
values Ai^5 of each lines are Ai = 0.3, A2 = —0.2, A3 = —0.34, A4 = —0.25, A5 = —0.2. The 
broken line shows the bounce solution. 




Fig. 10. The value of the smallest eigenvalue A of D as a function of (f)(0)- 
(1) The proper valley contains the bounce solution (broken line in Fig. 9). This solution 
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has a thick wall, as expected. It has a negative eigenvalue A3, which causes the instability. 
Another solutions that construct the valley are not the classical solutions. We call these 
series of the thick- wall (including the bounce) to the thin- wall bubbles "valley bubbles". (2) 
In Fig. [TO, when 0(0) goes to zero, A goes to one, because the minimum eigenvalue in the 



false vacuum is equal to the mass, which is equal to one. The eigenvalue A3 of the bounce 
solution does not occupy any particular position in Fig. [TO. As 0(0) approaches to one, A 



approaches to zero. This is the region where we have large thin-wall bubbles, as can be seen 
in the As-line of Fig. 9. (There is another zero point near 0(0) ~ 0.5, but this point does not 
play any special role.) The interior of these bubbles is the true vacuum, = 1. The latter 
property is especially notable: Even though the bounce solution is a thick-wall bubble, we 
find that the valley contains large, thin-wall, clean bubbles in the outskirts. 

These thin- wall bubbles can be analyzed by extending the original Coleman's argument: 
The solution of the proper valley equation Eq.(4-40) and Eq.(4-41) extremizes the action 
Spy + Sp. Using the valley equation Eq.(3-26), we rewrite it as follows for negative A; 

S P y + Sp = S + S x , S x = ^-Jd*xF 2 . (4-42) 

Now consider a fictitious particle in a one-dimensional space at time p. In the equation 
Eq.(4-40), the auxiliary field —F(p) acts as an "external force" for this particle. If the initial 
value 0(0) is sufficiently close to 1 and F is sufficiently small, remains close to 1 for a 
long time, until the friction term dies away. When it finally rolls down the hill, it does so 
under the external force —F(p). If —F(p) is just right, the particle approaches to the top of 
the lesser hill = asymptotically. Note that there is a major difference with Coleman's 
argument here: When the roll-down occurs, it does so under the external force, which is the 
sole source of the energy reduction, while in the bounce solution the timing of the roll-down 
has to be such that the friction term is just right to take care of the extra energy. [This does 
not prove the existence of the thin-wall solution, but Fig. 9 shows that this in fact happens.] 
Since —F(p) is the stopping force, it deviates from zero only at the wall. Therefore the 
second term in Eq. (|4-42j ) contributes positively only at the wall. The actions S and S\ can 
be approximately written in terms of the radius of the bubble R as the following; 

7T 2 

S = -ey# 4 + S 1 2n 2 R 3 , (4-43) 
Sx = ^vr 2 ^, (4-44) 

where Si and Wp are the numbers of 0(1). Taking the derivative of S + S> in the above 
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with respect to R, we find the radius of the solution of the proper valley equation to be 

Rrv = -(s 1 + ^-). (445) 



e 



2 A 



Therefore, even when e is not small enough to guarantee the large radius, |A| can be small 
enough to do so. This is what is causing the thin-wall bubble to be a solution of the proper 
valley equation. 

The shape of the large bubble and its thin wall can be examined in detail by looking 
at the 0+1 dimensional model, since the friction term can be neglected. We have thus also 
analyzed the 0+1 dimensional model. (See the previous subsection.) We have obtained the 
numerical value Wp ~ 0.2104 from this analysis. The radius seen in Fig. 9 for A45 is in 
agreement with Eq. Q4-45|) for this value of Wf- 



Numerical values of the action S are plotted in Fig. 11, where the horizontal coordinate 




Fig. 11. The action of the solutions of the proper valley equation as a function of the norm \<p\. 
The points with eigenvalues Ai^4 correspond to the lines in Fig. 9. 

is the "norm" of the solution, |0| = sj J d 4 x 4>(x) 2 , which for large R should be ~ ^Jn 2 /2R 2 . 
We see in this figure that the bounce solution denoted by its smallest eigenvalue, A3, is in 
fact at the maximum point of the valley. We have also compared the asymptotic expression 
of the action S in Eq.(4-44) with numerical value of S for large |0| and have confirmed that 
the action is in fact dominated by the volume and the surface terms. 

On the other side of the valley are small bubbles. In order to see their role, we have 
calculated the energy E of the mid-section of bubbles and plotted the result in Fig. |T2|. As 
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Fig. 12. The energy of the bubbles at the mid-section, where d(p/dr = 0. 

is well known, the bounce solution (A3) has E = 0, so that the energy conservation allows 
it to contribute to the tunneling from the false vacuum. We find that the bubbles smaller 
than the bounce have E > 0, while the larger ones have E < 0. This shows that the smaller 
bubble contribute to the tunneling from the states with higher energy than the false vacuum. 

We have numerically examined the valley that contains the bounce solution. We have 
obtained a series of deformed bounce solutions, valley bubbles, on the valley. These config- 
urations play the roles similar to instanton anti-instanton pairs with various distances. We 
have found that even when the bounce solution is a thick-wall bubble, the valley contains 
thin-wall large bubbles, whose interior is the true vacuum. Smaller bubbles are also identi- 
fied. We expect them to contribute to the decay of the higher energy states. This is quite an 
interesting possibility, which should be explored further. We note that we confined ourselves 
to spherical configurations in this letter. However, non-spherical configurations may also be 
interesting, as they could also contribute to the induced decay. 

4.3. Valley instanton in the scalar 4 field theory 

In the following two subsections, using the proper valley method, we examine the 
instanton-like configurations in the unstable scalar field theory and the gauge-Higgs sys- 
tem. il 1 In both systems, the instanton-like configuration has an important role in the quan- 
tum tunneling, but a simple scaling argument shows that no classical solution exists except 
a zero size one. 
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To construct finite size instanton-like configurations, it was often used the so-called the 
constrained instanton method, 0-* which was explained briefly in section 2 and 3. In this 
method, one introduces a constraint to define a sub-functional space of finite radius config- 
urations. The field equation is solved in this subspace and the finite-radius configuration 
similar to instantons is defined. After doing the one-loop (and possibly higher order) integral 
in this subspace, the constraint parameter is integrated over. The configuration constructed 
this way is called "constrained instanton" . One problem about this method is that its validity 
depends on the choice of the constraint: Since in practice one does the Gaussian integration 
around the solution under the constraint, the degree of approximation depends on the way 
constraint is introduced. Unfortunately, no known criterion guarantees the effectiveness of 
the approximation. 

This situation could be remedied once one realizes that what we have near the point- 
like (true) instanton is the valley.© That is, although the zero-size instanton may be the 
dominating configurations, it is expected to be followed by a series of configurations that 
makes the valley of the action. Therefore, instead of trying to cover all the neighborhoods 
of the zero-radius instanton, one may cover the valley region, which is expected to dominate 
the path-integral. The trajectory along the valley bottom should correspond to the scaling 
parameter, or the radius parameter of the instanton. As such, the finite-size instanton can 
be defined as configurations along the valley trajectory. This is similar to a calculation in the 
electroweak theory in which one evaluates the contribution of the instanton-anti-instanton 
valley. @>@)-@) Thus treating a single instanton as a configuration on the valley provides 
a means of unifying the approximation schemes. These configurations are named "valley 
instanton" . 

One convenient way to define the valley trajectory is to use the proper valley method. 
Since the constrained instanton has been used extensively in the existing literature, our main 
purpose here is to establish the existence of the valley instanton on a firm basis and compare 
its properties with that of the constrained instanton. This should serve as a starting point 
for the reanalysis of the existing theories and results under the new light. 

First, in this subsection, we consider the scalar 4 theory in the 4 dimensional space-time 
with Euclidean Lagrangian, 



£ = 1 

g 2 



(4-46) 



The negative sign for the 4 term is useful for estimating the asymptotic behavior of the 
perturbative expansion in the theory with a positive 4 term. 0^0* >H>>il)>ii3) Namely, the 
instanton-like configurations in the negative 4 theory reveal the large order behavior of the 
positive 4 theory. This model allows finite-size instantons only for /x = 0; otherwise, a 
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scaling argument shows that the action of any finite-size configuration can be reduced by 
reducing its size. As is shown in the following, when the size of the configuration is small, 
the valley instanton can be constructed analytically. The differences from the constrained 
instanton are rather small in this case. For a configuration with a large radius, we construct 
both instanton numerically. It will be found that remarkable differences between them 
appear in this case. While the action of the valley instanton increases monotonously and 
remains positive, the one of the constrained instanton becomes negative when the radius 
becomes large. 

4.3.1. Valley instanton 

For this model, the proper valley equation introduced in the previous section is 

-^0 + /i 2 0-^0 2 = F, 



(4-47) 

(-<9A + /i 2 - -0 2 ) F = XF. 



Now we introduce the scale parameter p defined by 0(0) = 4v / 3/p, in order to fix the radius 
of the instanton solution to be unity. We rescale fields and variables as in the following; 



Jl X = p 2 u, 0(x) = ^, F(x) = £f(r). (4-48) 



P P P 

The proper valley equation (4-47) is given by the following under this rescaling; 

+ {ppfh - i/* 3 = (p/i) 2 /, (4-49) 

+ (ppff - h 2 f = {ppfvf. (4-50) 

This system has an instanton solution in the massless limit, pp — > 0.0^'© In this limit, 
(4-49) reduces to the equation of motion and (4-50) the equation for the zero-mode fluctua- 
tion, /, around the instanton solution. The solution is the following; 

• ^ A = ^_ (4 . 51) 



1 


d 






dr 


1 


d 






dr 



u 1 + r 2 ' JV \l+r 2 (l + r 2 )\ 

where C is an arbitrary constant. Note that solution f is obtained from d(p (x)/dp, O — 
ho/ p. 

Let us construct the valley instanton in the scalar 4 theory analytically. When pp is 
very small but finite, the valley instanton is expected to have pp corrections to (4-51). On 
the other hand, at large distance from the core region of the valley instanton, since the term 
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of 0(h 2 ) is negligible, the valley equation can be linearized. This linearized equation can 
be solved easily. By matching the solution near the core and the solution in the asymp- 
totic region in the overlapping intermediate region, we can construct approximate solution 
analytically. We will carry out this procedure in the following. 

In the asymptotic region, (pp) 2 h 2 , the linearized valley equation is the following; 



1 d 

r 3 dr 

1 d 

r 3 dr 



,dh 



r 



r 



The solution of these equations is 



dr 
dr) 

f 



! + (pp) 2 h = (pp) 2 f, 
+ (ppff = {ppfvf- 



(4-52) 



(4-53) 



h(r) = C 1 G p ,(r) + J -, f(r) = C 2 G p ^(r), 
where C\ and C 2 are arbitrary functions of pp. The function G m (r) is 

m{r) = (27r) 2 r ' 

where K\ is a modified Bessel function. The functions / and h decay exponentially at large 
r. In the region of r <C (pp)" 1 , r <C {pp\/l — v) -1 , / and h can be expanded in series as the 
following; 



(4-54) 



h 



+ 



/ = 



Ci 
(2vr) 2 

C 2 
(2tt)V 

C 2 

(2tt) 2 



+ o(P/^) 2 ln(p/xrc) + 



^2 + 2^ 2 ^ ~ ^ ln(p/iVl - vrc) + 
"2 + ^(P^) 2 ( X - ^ HPP^l - vrc) + • 



(4-55) 



In the above, c = e 7_1//2 /2, where 7 is the Euler's constant. 

Near the origin, we expect that the valley instanton is similar to the ordinary instanton. 
It is convenient to define h and / as the following; 



h = h + (ppfh, f = f + (ppff, 



(4-56) 



where C in / is the function of pp and is decided in the following. The "core region" is 
defined as 

h > (pp) 2 h, f > {ppff. (4-57) 
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The valley equation for perturbation field h, f becomes 

-hi ( r3 §) " \ hlh = /o " ,! °- 

) . (4-58) 
The left-hand side of (4-58) has a zero mode, <p. It satisfies the equation, 



and is given by the following; 

Ay/3 8^ 
1 + r 2 (1 + r 



We multiply r 3 (p to both sides of (4-58), and integrate them from to r. The existence of 
zero mode ip makes it possible to integrate the left-hand side of (4-58). As a result of the 
integration by parts, only the surface terms remain and we obtain, 

dh 3 dip 



r 3 ™ + = f dj j r * r /o _ (4 . 61) 

dr dr Jo 



_ r^^L + r 3^/ = f dr y3 [(* - l)/o + • (4-62) 
dr dr Jo 1 J 

First, using (4-51) and ( |4-60|) , we can find that the right-hand side of (4-61) is proportional 
to lnr at r ^> 1. Thus (4-61) becomes 

dh 

r— + 2h = 4:V3(l-C)lnr. (4-63) 
dr 

This equation can be solved easily at r ^> 1. Using the solution of (4-63), we can find the 
right-hand side of (4-61) also proportional to lnr as the following; 

r < ^-+ 2f = iy/SC(l -u) lnr. (4-64) 
dr 

Finally, we obtain h and / at r 3> 1, 

h = 2y/3(l -C)lnr + / = (1 - u)2y/3 Clnr + • • • . (4-65) 
For these solutions to meet to (4-55), the parameters need to be the following; 

d = 0, C 2 = 4v / 3(2vr) 2 , u = 1, C = 1, (4-66) 
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as pp = 0. Now we have obtained the solution of the proper valley equation; 
4x/3 



h(r) 



1 + r 2 ' 
4V3 



4\/3(27r 2 )G 



if r <C (pp) 



-1/2. 



o((pp) : 



if (pp)~ 1/2 < r < (pp) -1 ; 
if (pp)~ 1/2 < r. 



(4-67) 



Let us discuss the consistency of our analysis. In the construction of the analytical 
solution, especially in the argument of the matching of the core and asymptotic region 
solution, we have implicitly assumed that there exists an overlapping region where both 
(4-52) and (4-58) are valid. Using the solution ([4- 671) , it is found that (4-52) is valid in the 
region of r ^> (pp) -1 / 2 , and (4-58) is valid in the region of r <C (pp) -1 . Therefore in the above 
analysis we have limited our calculation in the overlapping region (pp) -1 / 2 <r <C (pp) -1 . 

We calculate the action of the valley instanton using the above solution. Rewriting the 
action in terms of h(r), we find 



S 



7l~ 

g 2 Jo 



drr i 



dr 



+ (pp) 2 /, 2 - ^ 



(4-68) 



Substituting the analytic solution for S, we obtain 

16tt 2 



S 



oUppY 



(4-69) 



The leading contribution term comes from the ordinary instanton solution, and the correction 
term comes from the distortion of the instanton solution. 

4.3.2. Constrained instanton 

Here, we consider the constrained instanton in the scalar (p 4 theory, following the con- 
struction in Ref. |10]). We require the constraint in the path- integral. The field equation 
under the constraint is 



SS 50 n 
— + 0-— = 0, 



(4-70) 



where a is a Lagrange multiplier. The functional O had to satisfy the certain scaling prop- 
erties that guarantee the existence of the solution. 0) We choose it as follows; 



O 



d x 



(4-71) 
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This choice is one of the simplest for constructing the constrained instanton in the scalar 4 
theory. Again, adopting the rescaling (|4-48|) , the equation of motion under this constraint 
becomes 



(4-72) 



where we rescale the parameter a as a = (pfi) 2 a. The solution of this equation can be 
constructed in a manner similar to the previous subsection. To carry out the perturbation 
calculation in the core region, we replace the field variable as h = h + (pfi) 2 h, where 
ho ^> (pp) 2 h. The field equation of h becomes 



r s dr\dr 2 



(4-73) 



We multiply this equation by the zero mode ( 4-6U| ) and integrate this from to r. Then we 
obtain 

o dh ndip f f T 
- r <p— + r A -^-h = - 
dr dr Jo 

The solution of this equation in the region where r 3> 1 is 



dr'r' 3 ip ho + &h\ 



(4-74) 



h = 2V3\nr- 1 -^a + .. 



(4-75) 



In the asymptotic region, we consider the field variable and the parameter as h 2 <C (pp) 2 , 
ah 4 <C 1. Under this condition, the field equation of the asymptotic region becomes 



1 d ( 3 dh" 
r 3 dr \ dr 



+ (pfifh = 0. 



(4-76) 



The solution is 



h{r)=C x G Mi {r), 



(4-77) 



where C\ is arbitrary constant. In the region of r (pp) , this solution can be expanded 
as the following; 



h 



(27T) 



"a + ni.Pl 1 ) 2 ln(pprc) + 



(4-78) 



Matching this solution and the core region solution (|4-75 ), parameters are determined as the 
following; 



d = 4^(2^)2, a 



96 



ln(pp). 



(4-79) 
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To summarize, the analytical solution of the constrained instanton we have obtained is 
the following, 

4x/3 



h(r) 



1 + r 2 



if r <C (pp)" 



-1/2. 



2\/3(pp) 2 ln(pprc) + o ((pp) 2 ) , if (pp)~ 1/2 < r < (p/x)" 
4V3(27r 2 )C pM (r), if (pp)- 1 /* « r . 



(4-80) 



The action of the constrained instanton is given by the following; 



S = ^T" - 967r2 (^) 2 ln (P^) + ( W 2 ) • ( 4 - 81 ) 

This differs from the action of the valley instanton at the next-to- leading order. This correc- 
tion term shows that the constrained instanton is more distorted from the ordinary instanton 
than the valley instanton. 

4.3.3. Numerical analysis 

In this subsection, we calculate the valley equation (4-49), (4-50) and the constrained 
equation (4-72) numerically. Then we compare the valley and the constrained instanton. 

Each of the equations is the second order differential equation, so we require two boundary 
conditions for each field variable to decide the solution. We require all the field variables 
are regular at the origin. The finiteness of the action requires h, f — > faster than 1/r 2 at 
infinity. In solving (4-49) and (4-50), we adjust the parameter v and f(0) so that h, f — > 
at infinity for the fixed pp. In the similar way, in case of the constrained instanton, the 
parameter a is determined so that h — > at infinity. 

Numerical solutions of the valley equation near the origin are plotted in Fig. 13 (a) for 
pp = 0.1, 1.0. The solid line shows the instanton solutions (4-51), which corresponds to 
pp = 0. The numerical solutions of the constrained instanton are also plotted in Fig. 13 (b) 
for pp = 0.1, 0.5, 1.0. Both the valley and the constrained solution for pp = 0.1 agree with 
the analytical result. As pp becomes large, both solutions are deformed from the original 
instanton solution. We find that the distortion of the constrained instanton is much larger 
than that of the valley instanton. This also agrees with the analytic result. In the analytical 
solution ( |4-80|) , the correction term 2v^3(pp) 2 ln(pprc) contributes to this distortion. On the 



other hand, the correction term of the valley instanton (|4-80| ) is o((pp) 2 ), which is smaller 
than the previous one. In addition, we find that the exponentially damping behavior of the 
analytical solution in the asymptotic region, where r ^> (pp) -1 agrees with the result of 
numerical analysis. 
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(a) 



(b) 




Fig. 13. (a) Shapes of the numerical solution of the valley instanton, h(r), for pp =0, 0.1, and 1.0 
near the origin. The solid line denotes the original instanton, ho- (b) Shapes of the constrained 
instanton for pp =0, 0.1, 0.5, and 1.0. 
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Fig. 14. (a) The action S (in units of <? 2 /167r 2 ) of the numerical solution of the valley equation as 
a function of the parameter pp. (b) The action S of the constrained instanton. 



The values of the action of the valley instanton for pp =0.0001 ~ 1.0 are plotted in 
Fig. 14 (a). If pp is very small, this result is consistent with ( |4-69| ). The values of the action 
of the constrained instanton for pp = 0.0001 ~ 0.8 are plotted in Fig. 14 (b). This figure 
shows that the behavior of the action is similar to that of the valley instanton when p\i is 
very small. When pp becomes large, the behavior of the action is different from the valley 
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instanton case. 



4.4. Valley instanton in the gauge-Higgs system 

In this subsection, we consider the SU(2) gauge theory with one scalar Higgs doublet, 
which has the following action S — S g + S^, 

S 9 = ^J d * x K F v»> F v» ( 4 ' 82 ) 
S h = \j d*x{(D,H? (D,H) + 1 (&H - v 2 ) 2 } , (4-83) 

where = d^A u — d^A^ — i L4 M , A v \ and = — iA^. The masses of the gauge boson 
and the Higgs boson are given by, 

1 

mw = V 2A^' m " = V2 V ' ^ 4 ' 84 ^ 

When v is not zero, this theory is known not to have any finite-size instanton solutions, 
in spite of its importance. In the following, we will construct instanton-like configurations 
for this theory, relevant for tunneling phenomena, including the baryon and lepton num- 
ber violation processes. As well as the unstable scalar theory, the valley instanton can be 
constructed analytically when the size of the configuration is small. For a configuration 
with a large radius, numerically constructions are done. There is differences between the 
valley instanton and the constrained instanton in this case. Therefore, when the large size 
configuration dominate the path integral, we must use the valley instanton instead of the 
constrained instanton. 

4.4.1. Valley instanton 

The valley equation for this system is given by, 

* 2 -F A + J^L;F H + JXt-F h ^ = X e F" 



M 5A^ " 6H' 

where the integration over the space-time is implicitly assumed. The valley is parameter- 
ized by the eigenvalue A e that is identified with the zero mode corresponding to the scale 
invariance in the massless limit, v — > 0. 
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To simplify the equation, we adopt the following ansatz; 



AJx) = ^^--2a(r), H(x) = v (1 - h(r)) rj, (4-86) 
x 2 

where rj is a constant isospinor, and a and h are real dimensionless functions of dimensionless 
variable r, which is defined by r = \fx 2 / p. The matrix a^ u is defined, according to the 
conventions of Ref. |6|), as a^ v = fj afiv a a /2. We have introduced the scaling parameter p so 
that we adjust the radius of valley instanton as we will see later in subsection 4.4.3. The 
tensor structure in (4-86) is the same as that of the instanton in the singular gauge. 

Inserting this ansatz to (4-85), the structure of and F H ^ is determined as the following; 

Ff{x) = • ^fr(r), F B \x) = ~f h (r) V . (4-87) 

X A A 

By using this ansatz, (4-86) and (4-87), the valley equation (4-85) is reduced to the following; 

-\tt ( r t) + v^ a - 1)(2a - 1} + fx (H2fl(1 - hf = x (/w)2/a ' (4 ' 88) 
( r3 S) + ^ - 1)fl2 + \ {pv)2h{h -^ h -v = ( 4 - 89 ) 

r dr \ dr J r z IX 

+ 9 l{ pv fa{h - l)f h = ^(pv) 2 uf a , (4-90) 
1 d ( r ^) + ^f + h pv y {3h 2_Q h + 2)fh 



r 



3 dr \ dr I r 2 4 



Ah-l)f a = {pvfvf\ (4-91) 



where v is defined as A e = v 2 u/\. 

In the massless limit, pv —>■ 0, (4-88) and (4-89) reduce to the equation of motion and 
(4-90) and (4-91) to the equation for the zero-mode fluctuation around the instanton solution. 
The solution of this set of equations is the following; 

i / 2 \ V 2 

a o — z — ; — nj h = 1 



1+r2 vi+^v ^ (4 . 92) 

fa _ 2Cr 2 h = Cr 

Jo (1 + r 2 ) 2 ' Jo (1 + r 2 ) 3 / 2 ' 

where C is an arbitrary function of pv. Note that ao is an instanton solution in the singular 
gauge and h is a Higgs configuration in the instanton background.!) We have adjusted the 
scaling parameter p so that the radius of the instanton solution is unity. The mode solutions 
/o and /o are obtained from dao/dp and dho/dp, respectively. 
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Now we will construct the valley instanton analytically. When pv = 0, it is given by 
the ordinary instanton configuration ao, ho, /q and f{f. When pv is small but not zero, it 
is expected that small pv corrections appear in the solution. On the other hand, at large 
distance from the core of the valley instanton, this solution is expected to decay exponentially, 
because gauge boson and Higgs boson are massive. Therefore, the solution is similar to 
the instanton near the origin and decays exponentially in the asymptotic region. In the 
following, we will solve the valley equation in both regions and analyze the connection in 
the intermediate region. In this manner we will find the solution. 

In the asymptotic region, a, h, f a and f h become small and the valley equation can be 
linearized; 



1 d 

r dr 



da 
dr 



4 

H — t;CL 



9_ 
2A 



(pvy a = »-(pvyr, 



-A{r"fV\{pv?h = {pvy f \ 



r 3 dr 
1 d I 
r dr \ 
]_d_ 

r 3 dr 



dt\ 
dr J 

dr 



+ 



1 



A 



f ^(pv) 2 f h = (pv) 2 uf h . 



(4-93) 
(4-94) 
(4-95) 
(4-96) 



The solution of this set of equations is 



d 



< r ) = Cir—G fmi J 



h{r) = 

f a {r) 
f h (r) 



C 2 G pmH (r) + -f h (r) 



C 3 r^-G 
dr 



1 

(r), 



C 4 G p , H (r), 



(4-97) 

(4-98) 

(4-99) 
(4-100) 



where Ci are arbitrary functions of pv and p, w H are defined as p w H = m w H y/l — 2u. As 
was expected above, these solutions decay exponentially at infinity and when r <C (pv)^ 1 
they have the series expansions; 



a(r) = 
h{r) = 

f a {r) - 



C l 

(27T) 2 
(27T) 2 



C 3 

(27T) 5 



2 !/ ^2 



+ 



c 3 



z/(2tt) ; 



2 !/ 
~2 + ~ 2 ^ 



+ 



1 



+ n(P m H?HP m H rC ) + 



c 4 



z/(2vr) 



^2 + \ (PPH) 2 HPP H rC ) + 



~l2 + «(PPw) 2 + 



(4-101) 

(4-102) 
(4-103) 
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f\r) 



c. 



(27T) 



(4-104) 



c being a numerical constant e 7 ~ ' /2, where 7 is the Euler's constant. 

Near the origin, we expect that the valley instanton is similar to the ordinary instanton. 
Then the following replacement of the field variables is convenient; a = a,Q + (pv) 2 a, h = 
h + (pv) 2 h, f a = fS + (pv) 2 f a , f h = / h + (pv) 2 f h . If we assume a > (pv) 2 a, h > (pv) 2 h, 
/o 3> (pv) 2 f a and /q ^ (pv) 2 f h , the valley equation becomes 



'0 • 



(4-105) 



M ( r3 f ) + ^ + r^° " lWa + \ ho{ho ~ 1){ho ~ 2) = fl (4406) 



1 d ( df aS 
r dr \ dr 



24 



+ -(6a 2 - 6a + 1)/- + — (2a - l)/ a a 



1 d 

r 3 dr 



dr 



+!L(h - l) 2 / a + ja Q (h - 1)/* = L vfSj (4-107) 
+ ^ a o/ ft + ^0/0* + |(3^g - + 2)/* 



r> 

+ -a (/i - l)/ a + ^0 - 1)/qS + ^«o/o ^ = Wo • (4-108) 

To solve this equation, we introduce solutions of the following equations; 

4 



l_d_ ( dl Pa\ 
r dr \ dr J 

3 d(f h 



+ ^(6a - 6a + l)</?a = 0, 



r 



(4-109) 



1 d 

r 3 dr 



r 



dr 



+ -^a </? ft = 0. 



They are given as, 



1/2 



(fa 



1 + r 2 



(4-110) 



(l + r 2 ) 2! 

Using these solutions, we will integrate the valley equation. We multiply (4-105) and (4-107) 
by np a , and multiply (4-106) and (4-108) by r 3 (fh then integrate them from to r. Integrating 
by parts and using (|4-109|) , we obtain 



da dtp a A g 2 r , , 
-p a r— + -j—ra = — dr r ip a 
dr dr A Jo 



fS - \ao(ho ~ I) 2 



(4-111) 



,dh 
dr 



dr 



dr'r' 3 (p h 



fo ~ \h (h - l)(h - 2) - ^(h - l)a a 



(4-112) 
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df a dcp a , 

-Var-j- + -j-rf a 



dr dr 

rr 

= / dr'r'(p a 
Jo 



9 
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^-ufS - - 2 {2a - l)fSa - ^(h - lYfS - »-a (h - 1)/* 



2A 



A 



(4-113) 



: ''" J lir'' (,r ' h 



J 

Jo 



•dj^ 
dr dr 

dr'r' 3 ifih 



vft ~ 4«o/o« - \{Zhl - 6d + 2)/ ft 



~a (h - l)/ a - 4(^o " l)/ °a " ^offh 



(4-114) 



First we will find a. The right-hand side of (4-111) is proportional to (C — 1/4) and when 
r goes to infinity this approaches a constant. At r 3> 1, (4-111) becomes 

Ida 2 



r dr 



Then at r ^> 1, d(r) is proportional to (C — l/4)r 2 and a(r) becomes 



1 



(pvfg 2 f c 1 
12A 



r 2 + 



(4-115) 



(4-116) 



To match this with (4-101), it must be hold that C = 1/4 when pv = 0. When C = 1/4, 
the right-hand sides of (4-111) vanishes and d satisfy —ip a dd/dr + ddip a /dr = 0. Hence d is 
d = -D <^a, where D is a constant. Identifying a + (pv) 2 d with (4-101) again at r ^> 1, we 
find that C\ + C 3 /f = — 27T 2 and C 3 = — 7r 2 at pv = 0. In the same manner, h, / a and f h 
are obtained. At r 3> 1, we find 



h = const. + 

to ^ ^ 

f a 



f h 
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r | r 2 - ^(1 -2v) + 



16A 



(4-117) 



1 - 2u) lnr + • 



Here const, is a constant of integration. Comparing (4-102)-(4-104) with them, we find that 
C 2 + C A /u = 2tt 2 , C 4 = vr 2 and v = 1/4 at pv = 0. 

Now we have obtained the solution of the proper valley equation. Near the origin of the 
valley instanton, r <C (pm WH )~ 1 ^ 2 , it is given by, 



a(r) = 



„2 ' 



2(1 + r 2 ) 2 ' 



h(r) = 1 
f h (r) = 



r 



1/2 



1 + r ' 2 , 

r 



(4-118) 



4(l +r 2)3/2' 
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where we ignore the correction terms that go to zero as pv — > 0, since they are too small 
comparing with the leading terms. As r becomes larger, the leading terms are getting smaller 
and so the correction terms become more important; 



1 1 (nv) 2 

a(r) = 72 + o( (p t ,) 2 ) , Mr) = _ - i^L 1„ 2 + ■ 



(4-119) 



for (pm WH ) ^ 2 < r « (P m W n) Finally, far from the origin, r ^> (pm WH ) 1//2 , the 
solution is given by the following: 

Mr) = -27T 2 G pm (r) + -f h (r), 

V J pm H \ , \ J, (4-120) 

/ ft (r) = 7 r 2 G p ^(r), 

where u — 1/4 for pt> = 0. In ( [4-1 20 ), we ignore correction terms, since they are too small. 

Let us make a brief comment about the consistency of our analysis. Until now, we have 
implicitly assumed that there exists an overlapping region where both (4-93)- (4-96) and 
(4-105)-(4-108) are valid. Using the above solution, it is found that (4-93)-(4-96) are valid 
when r 3> (pm w H )~ 1 / 2 and (4- 105)- (4- 108) are valid when r <C (pm WH )~ 1 . If pm WH is 
small enough, there exists the overlapping region (pTn WH )^ 1 ^ 2 <r<C (pm WH )~ l . Then our 
analysis is consistent. 

The action of the valley instanton can be calculated using the above solution. Rewriting 
the action in terms of a and h, we find 

Yin 2 r 00 dr \ ( da x 2 



S„ = ^(P«) 2 [ r 3 dr | (f) 3 + ^(h - 1)V + \(ftv?h 2 (h - 2)4 . (4.122) 
Substituting the above solution for S, we obtain 

S = ^ + - ^( P v) 4 Hpv) + 0((pv)% (4-123) 

The leading contribution 87r 2 /g 2 comes from S g for ao, which is the action of the instanton, 
and the next-to- leading and the third contributions come from Sh for ao and ho. 
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4.4.2. Constrained instanton 

Here, we consider the constrained instanton. According to Affleck's analysis, 0* the 
constrained instanton satisfies the following equation: 

SS 50 A 

'a-/ = 0, (4-124) 



5A„ 5A^ 

where a is a Lagrange multiplier and depends on the constraint. Both Oa and Oh are 
functionals of and H respectively that give a solution of the constrained equation, as the 
scalar theory in the subsection 4.3.2. Here we adopt the ansatz (4-86) again. By the similar 
analysis as the valley instanton, it turns out that the behavior of the constrained instanton 
is the following. Near the origin of the instanton, the solution is given by Oq, h Q in (4-92) 
as well as the valley instanton. In the region where (pm WH )~ 1 / 2 < r <C {P m w h)~ 1 ■> the 
solution is given by 

1 (pv) 2 q 2 1 (pv) 2 

a(r) = - - ffi-L. + . . • , h{r ) = — + ^- In(pvrc) + ■■■. (4-126) 

Let us compare ([4-1 19|) and ([4-126|) . The correction term of the valley instanton is smaller 
than one of the constrained instanton. Finally, for r ^> (pm w H )^ 1 ^ ! , the constrained instan- 
ton is given by 

a(r) = -^r^-G^r), h(r) = 27r 2 G pm ». (4-127) 
The action of the constrained instanton is the following: 

s = 5~ + T~ (/w)2 + 0((/w)4 ln(/w)) - (4 ' 128) 

The leading contribution 87r 2 /g 2 comes from S g for oq, which is the action of the instanton, 
and the next-to-leading contributions comes from Sh for ao and Hq. The difference from 
the valley instanton is that we cannot determine the term of 0((pv) A hi(pv)) by the current 
analysis. 

Now we choose a constraint and analyze the constrained instanton . We adopt the fol- 
lowing functionals for the constraint: Oa = ig 2 J d A xtiF^ v F up Fp^, Oh = 0. This constraint 
is one of the simplest for giving the constrained instanton. Then the constrained equation 
of motion is given by, 
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+6{pvf 



(2a +2a(a-l 



-a 2 (a-l) 2 (2a-l) 



d 2 a 



a(a — 1' 



da 



(4-129) 



0, 



■hi ( r3 f ) + ^ - 1)fl2 + i^ 2 ^ - - 2) = °« 



(4-130) 



where a = p A v 2 a. By solving these equations approximately, in fact, we obtain the behavior 
of the solution that we have given previously, and the multiplier is determined by a = 
5g 2 /A8\ as pv -> 0. 

4.4.3. Numerical analysis 

In this subsection, we solve the valley equation (4-88)-(4-91) and the constrained equation 
(4-129), (4-130) numerically, and compare the valley instanton and constrained instanton. 

We need a careful discussion for solving the valley equation (4-88)-(4-91): Since the 
solution must be regular at the origin, we assume the following expansions for r <C 1; 



a{r) = J2 a(n)r n , h(r) = E V) r ™> 

n=0 n=0 

oo oo 

f a (r) = E /(V"' /V) = E f(nf n - 

n=0 n=0 

Inserting (4-131) to (4-88)-(4-91), we obtain 

O(o) = 1, h( ) = 1, / ( a 0) = 0, / ( o) = 0, 
a (1) = 0,/ ( a 1) = 0, 

^(2) = 0,/f 2) = 0. 



(4-131) 



(4-132) 



The coefficients a(2), /(2) and are not determined and remain as free parameters. 
The higher-order coefficients (n > 3) are determined in terms of these parameters. Four 
free parameters are determined by boundary conditions at infinity. The finiteness of action 
requires a, h — > faster than 1/r 2 at infinity. This condition also requires f a , fh —> 0. 

We have introduced p as a free scale parameter. We adjust this parameter p so that 
a(2) = —2 to make the radius of the valley instanton unity. As a result we have four 
parameters h^, ffy, and pv for a given v. These four parameters are determined so 
that a, h, f a , and f h — * at infinity. In the case of the constrained instanton, the two 
parameters h(\) and a are determined under a (2) = —2 so that a, h — > at infinity. 
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A numerical solution of the valley equation near the origin is plotted in Fig. 15 for 
pv = 0.1, 1.0 at X/g 2 = 1, when m w = m H . We plot the instanton solution (4-92) by 
the solid line, which corresponds to the valley instanton for pv = 0. This behavior of the 
numerical solution for pv = 0.1 agrees with the result of the previous subsection, (4-118). 
Moreover, even when pv = 1.0, the numerical solution is quite similar to the instanton 
solution. We find that the behavior of f a (r) and f h (r) also agrees with the analytical result 
(4-118) as well as a(r) and h(r). We also find that the numerical solution in the asymptotic 
region where r (pt>)~ 1//2 , is damping exponentially and agrees with the analytical result 

(EH). 




Fig. 15. Shapes of the numerical solution of Fig. 16. Shapes of the numerical solution of 

the valley equation, a(r) and h(r) for pv = the constrained instanton, a(r) and h(r) 

0.1, 1.0 near the origin. The solid lines de- for pv = 0.1, 0.5, 1.0 near the origin. The 

note the original instanton solution, ao and solid lines denote the original instanton so- 

/i . lution, ao and /iq. 



On the other hand, a solution of the constrained equation is plotted for pv = 0.1, 0.5, 1.0 
at X/g 2 = 1 in Fig. 16. The numerical solution for pv = 0.1 also agrees with the analyt- 
ical result. As pv is larger, both the valley instanton and the constrained instanton are 
more deformed from the original instanton solution. Nevertheless, the correction of the con- 
strained instanton from the original instanton is much larger than one of the valley instanton, 
especially when pv>0.5. 

The values of the action of the valley instanton for pv = 0.001 ~ 1.0 are plotted in 
Fig. 17: Fig. 17(a) depicts the behavior of the total action S, while the contribution from 
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g 2 S/87t 2 
1.02 

1.01 



1 
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pv 




pv 
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Fig. 17. (a) The action S (in unit of g 2 /8ir 2 ) 
of the numerical solution of the valley 
equation, at X/g 2 = 1, clS EL function of 
the parameter pv. (b) The contribution 
from the Higgs sector, Sh- The solid line 
shows the behavior of the analytical result 
that g 2 S h /87r 2 = (/w) 2 /4-(H 4 ln(H/32. 
(c) The contribution from the gauge sector, 
S a . 



(a) 



(b) 



(c) 



g 2 S/Sn 2 
1.02 

1.01 
1 



pv 




pv 



pv 

0.5 1 

Fig. 18. (a) The action S (in units g 2 /8ir 2 ) of 
the numerical solution of the constrained 
instanton, at X/g 2 = 1, as a function of the 
parameter pv. (b) The contribution from 
the Higgs sector, Sh- The solid line shows 
the behavior of the analytical result that 
g 2 S h /8ir 2 = (pv) 2 /4. (c) The contribution 
from the gauge sector, S g . 



the gauge part S g and from the Higgs part Sh are in (b) and (c) respectively. The solid line 
shows the behavior of Sh of the analytical result ( |4-123|) . From Fig. 17 (c), it turns that 
S g is almost independent of pv. This figure shows that our numerical solutions are quite 
consistent with ( [4-123]) , even when pv = 1.0. 

The values of the action of the constrained instanton for pv = 0.001 ~ 1.0 are plotted 
in Fig. 18. When pv >0.5, we can no longer use the analytical result ( |4T28| ). This is 
attributed to the large deformation from the original instanton. We notice that the action of 
the constrained instanton is smaller than that of the valley instanton for the same value of pv. 
This result is natural at the point that the parameter pv corresponds to the scale parameter. 
This point was elaborated upon in section 3, using the toy model (3-23). To repeat, the fact 
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that the constrained instanton has the smaller action than the valley instanton for the same 
value of pv does not mean that the constrained instanton gives more important contribution 
to the functional integral than the valley instanton. 



We summarize all the numerical data in Table [TJ and III, for the valley instanton and for 
the constrained instanton, respectively. 



pv 


V 


h (1) 


/(2) 


rh 
1(1) 


9 2 S/8n* 


0.001 


0.2500 


-1.000 


0.5001 


0.2500 


1.000000 


0.005 


0.2500 


-1.000 


0.5001 


0.2500 


1.000006 


0.01 


0.2501 


-1.000 


0.5001 


0.2500 


1.000025 


0.02 


0.2501 


-1.000 


0.5003 


0.2501 


1.000100 


0.05 


0.2503 


-1.000 


0.5014 


0.2505 


1.000625 


0.1 


0.2509 


-1.000 


0.5044 


0.2515 


1.002503 


0.5 


0.2537 


-1.007 


0.5478 


0.2669 


1.062550 


1.0 


0.2521 


-1.021 


0.6226 


0.2926 


1.244646 



Table II. The numerical data of the valley instanton in the gauge-Higgs system. 
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pv 


a 


h(i) 


Pov 


g 2 S/87r 2 


0.001 


0.1042 


-1.000 


0.000100 


1.000000 


0.005 


0.1042 


-1.000 


0.005000 


1.000006 


0.01 


0.1041 


-1.000 


0.010000 


1.000025 


0.02 


0.1040 


-1.000 


0.019998 


1.000100 


0.05 


0.1035 


-1.002 


0.049972 


1.000623 


0.1 


0.1019 


-1.005 


0.099779 


1.002471 


0.5 


0.0798 


-1.070 


0.481738 


1.053903 


0.9 


0.0637 


-1.144 


0.833000 


1.152691 


1.0 


0.0609 


-1.162 


0.917406 


1.183336 



Table III. The numerical data of the constrained instanton in the gauge-Higgs sys- 
tem. 



In Table III, we show the parameter po used as "p" in Ref. [10]). The parameter po is 
defined by, 



O = c Po 



4-d 



(4-133) 



where d is a dimension of the operator and Co is a conveniently chosen constant. For Oa = 
ig 2 J d 4 xtYF fJiU F l/ pFp^ and Oh = 0, we choose that Co = 967r 2 /5g 2 so that po = p as pv — > 0. 
We find that po nearly agrees with p. 



4.5. Valley calculation of the B and L violation 

In this subsection, we calculate the cross section for the process of the type q+q — ► (3n/ — 
2)q+nfl+n w W+nf l H, using the valley instanton. 0) As explained in section 2, Ringwald and 
Espinosa studied this process by using the constrained instanton. They suggested that 
the amplitudes of the multi-boson process increase with the energy and reach lpb for a center 
of mass energy of O(10) TeV. Unfortunately, the high-energy behavior of the process with 
many bosons, which is most interesting, is out of the validity of their calculations. One of the 
problem that they used the analytically constructed configuration given in subsection 4.4.2. 
In the process with many bosons, the large radius instanton dominates the path-integral. 
As is shown in subsection 4.4.3, the large radius constrained instanton is not approximated 
well by the analytically constructed configuration. On the other hand, the valley instanton 
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can be approximated by the analytically constructed configuration even for the large radius 
one. (See subsection 4.4.3.) Thus the analytically constructed configuration can be used in 
the calculation. Although induced vertex by the valley instanton is similar to that given by 
Ref. 0) and [!]), the latter can not be justified in the calculation of the process with many 
bosons. 

Another problem of Ref. ^) is that the orientation effects of the instanton are not taken 
into account. The integration of the orientation in the Lorentz SU(2) R rotations is not trivial 
but very important especially for the process with many bosons. Moreover, in considering 
the multi-boson processes at high energies, masses of the gauge boson, Higgs particle and 
top quark can not be neglected. All the effects are incorporated in the following calculations. 

4.5.1. Valley instanton in the electroweak theory 

We consider the standard model with n/ generations of quarks and leptons. For simplicity, 
we set SU(3) C and U(l) y coupling constant to be zero. Thus we treat the simplified standard 
model that reduces to an SU(2) gauge-Higgs system with fermions. As for quarks, we have 
the left-handed doublets q l £ and the right-handed singlets q l R , where a = 1, 2, 3 is the color 
index and i = 1, • • • ,nf is the index of the generations. Similarly, we have the left-handed 
doublets l l L and the right-handed singlets l R for leptons. The bosonic part of the action is 
given by 

S 9 = 7T-o f d 4 x trF, u F, u , (4-134) 



2gK 

S h = \j d A x (D^H) + i (rfH - v 2 ) 2 } , (4-135) 

where F^ = d^A v — d u A^ — % [A^, A v ] and = — iA^. The masses of the gauge boson 
and the Higgs boson are 

mw/= V2A U ' m " = 72 V ' ^ 4 ' 136 ^ 
The fermionic part of the action is 

S q = J d A x(iq[a IM D ll q L - iu^a^UR - id R a^d^d R (4-137) 

+y u q[eH*u R - y d q ] L Hd R + h.c), 

Si = J d 4 x{il ] L a tl D tJl l Ij - ie^a^eR - yil ] L He R + h.c), (4-138) 

where e Q/3 = — e a p = ia 2 , <t m = (ex, i) and = (a, —i). The masses of fermions are given by 
m u,d,i = Vu,d,i v - The total action is given by 



S = S g + S h + y £[s l + y £S q ). (4-139) 
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As was shown in the previous subsection, the dominating configuration in the path- 
integral of the bosonic variables is the valley instanton. The valley instanton is given by the 
following; 

A^x) = %Ua^ ■ 2a(r), (4-140) 
H{x) = v(l - h(r))r], (4-141) 

where U is a Lorentz SU(2) R matrix, rj is a constant isospinor, and a and h are real dimen- 
sionless functions of dimensionless variable r, which is defined by r = \fx? / p. We choose r] r 
as (0,1). 

Near the origin of the valley instanton, r -C (pm WH )~ 1 ^ 2 , the valley instanton is given 

by 

a ^ = rh> W = 1 -(r^) 1/2 ' (4442) 

where we ignore the correction terms that go to zero as pv — > 0, since they are too small 

comparing with the leading terms. As r becomes larger, the leading terms are getting smaller 

and so the correction terms become more important; 

1 1 (nv) 2 

a(r) = - + o((pv) 2 ), h(r) = — -^-ln2 + ---, (4-143) 



r 



for (pm WH ) ^ 2 C r « (P m w.n) Finally, far from the origin, r ^> (pm WH ) 1//2 , the 
solution is given by the following: 



<r) = ^G pmw {r) + -nr) 
vr 2 „ , , 1 



2^ 2 dr v' ' " (4-144) 

Kr) = ~G, m (r) + -f h (r), 



2v 1 H " ' v' 

where v — 1/4 for pv = 0, and by the numerical calculation, we obtain v ~ 1/4 for pv < 1. 
In (|4-144|) , we ignore correction terms, since they are too small. The function G m (r) is 

G-M - (4.145) 

where i^i is a modified Bessel function. The functions / and h decay exponentially at large 
r. 

The fermionic zero mode around the valley instanton is given by; 



Tip- 

i m 



\x v aTUH (-e rs rf s VpMr) + t^e^d^r)) , (4-146) 



u Ra {x) = ^u R {r)Ul%, (4-147) 

Z7T p 

d Ra {x) = ±^d R (r)UVe fh rf', (4-148) 

Z7T p 
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where the Greek and the Roman letters denote indices of spinor and isospinor, respectively. 
For small pv, the solution is given approximately, by the following; 



and di and are obtained by replacing m u with mj. The above approximative behavior is 
correct even though pv is 0(1), since the valley instanton keep its shape unchanged as was 
seen in section 4.4.3. 

4.5.2. Calculation of the cross section 

To evaluate the multi-point Green function, we must introduce collective coordinates xq 
and U to restore the translational and Lorentz invariance as the valley instanton breaks 
these symmetries. In addition to them, the p integral must be brought into, because along 
the valley direction the action varies gently and the corresponding integration can not be 
approximated by the Gaussian integration. 

Since the valley instanton is not a classical solution, a slight complication arises to in- 
troduce the collective coordinates. For the classical solution, corresponding to the broken 
symmetries, there exist zero modes and the collective coordinates are introduced naturally 
to extract these zero modes from the functional integral. However, the existence of the zero 
modes is not guaranteed for the valley instanton. In Ref. 0), another collective coordinate 
method have been used. The point is that the collective coordinates for the translations 
and the Lorentz SU(2)ji rotations can be treated as transformation parameters of a kind of 
"gauge symmetries" . © The orbits of the "gauge symmetries" can be extracted by the usual 
Faddeev- Popov method and they becomes the collective coordinates for the translations and 
the Lorentz SU(2)r rotations. As well as the usual gauge symmetry, the translation and 
the Lorentz invariance of the Green function is automatically guaranteed. For the collective 
coordinate of the valley direction, we bring into it so as the quantum fluctuations around 
the valley instanton are restricted to be orthogonal to the valley direction. 0) From this 
prescription, liner terms of the quantum fluctuation drop naturally from the action while we 
do not treat the classical solution. 

Since the valley instanton is very similar to the original instanton, we approximate the 
determinant and Jacobian of the valley instanton by one of the original instanton!) . Then 




1/2 . 



1/2 . 



(4-149) 



(4-150) 
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the Fourier transformed Green function is given by, 

G(p,q,k) =e- 8 ^V 8 fduf dpe^ 2 ^ 2 / x p 2n f- 5 c( P fiY 43 ^/ 6 

x^) ■ ■ ■ A„( qi ) ■ ■ 4(h) ■ ■ ■ . (4-151) 

Here the Fourier transformations of the fermionic zero mode ip is given by 

H%)--HiK x )^ + - (4 - i52) 

where • • • includes the regular term in the limit: p 2 + m 2 — > 0, and \a is the constant spinor, 
(0, —1) for ip — u, v and (1,0) for ip — d,e. The Fourier transforms of the valley instanton, 
is given by 

A il = — 7T ip 2|2 H • (4-153) 

We define the shifted Higgs field (f> hy H = (v + (f))rj and the Fourier transform is given by, 

^"V^-- (4454) 



The integration of the valley parameter pv is given by 

jf 'w w/ vr = ^ (^V r (o, (4-155) 



where t = |n/ + n w + + The integration is dominated by the contribution around the 
saddle point; 

f \ 1 )!/2 

"Hs^-s*} (4 ' 156) 

Therefore, we can check the consistency of the approximation by this value of the saddle 
point. The saddle point pv s depends on the number of bosons monotonously. The valley 
instanton is quite similar to the original instanton even at pv = 1, and therefore the approx- 
imation is valid for n w + ~ 40, if we assume m H = m w and use g 2 ~ 0.42. On the other 
hand, the constrained instanton is deviated from the original instanton at pv = 0.5 and so 
the approximation breaks for n w + rih ~ 4. Therefore we can calculate the amplitude of the 
multi-boson process, using the valley instanton. Finally, we obtain the Green function; 

/ v(43-8n/)/6 

G(p, q, k) = HfV • 2- Uh/2 ^— m -6n f -2n w -n h+ 4 e ~8^/gl 

\ m w J 



x g Un f /3+n w +n h -29/6 J dUL (lJ. p> g> ^ (4457) 
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where 

c > = c . 2 -2n//3-25/6 7r -2n / /3-19/6j.^ 

= 1.26 x 10 6 • e 6 - 516n f ■ r(t), (4-158) 

and L(U ; p, q, k) is defined by, 



An, 



rrf-Y \ 1 ^ q iv Ua IHJ tf 2* 1 



We ignore the momentum dependence that goes to zero on mass-shell. Performing the LSZ 
procedure, we obtain the invariant amplitude T nj , nvj , nh . Summing up the polarization and 
charge of the gauge field, isospinor and spinor of the fermion field, we obtain the amplitude; 

/ „ \(43-8n/)/3 

E\ T |2 _ l2 - nu ,-n h I _V_ \ / 2 \-Qn f -2n w -n h +A -16ix 2 /g 2 

\J-n f ,n w ,n h \ — C Z ™ v'V/ 



4ny „, 

x ^2 ) 14n / /3 + n u ,+n h -29/6 / ffl J| ~Q (-)gW q » q ° t r(<J^U<T p *tf) . 

i=l j=l 

(4-160) 



The cross section is given by the following; 

a = d^i h / 1 (4^ n (2^32^. n (2^^ (27r)4 ^ (4) fen - 

X E l^.n^.nj 2 • , 1 (4-161) 

4y(pi -P2) 2 - m\m\ 

J2 / \(43-8n/)/3 

xe -16.W (^14^/3+^-29/6 . J^^) . ,-1, 

where s = pf n = (pi + p<i) 2 and the masses of fermions of the initial state are ignored. We 
denote the sum of the momentums of the final state by p out . All the information about phase 
space and group integration is encoded in the function l2n f ,n w ,n h ( s )- We define //, m , n (s) 

r 21 d 3 p- m d 3 q- n d 3 k[ 

h,m,n( S ) = / II ^yp~ II II ^"^(Pin ~ Pout) 

J i=3 j=l J=i ltj h 

21 m 

x2 m f dufltrthU) f[(-)9^q]tT(a UI/ Ua pa U^ (4-162) 

J i=l 7=1 



J i= 3 ^Pi 7=1 J=l ^fc, 



22 



/^Jb (ft 
d^ntaWOlI Hi^)} 

i=l 7 = 1 L 



2 2 
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To evaluate the function /i, m , n (s), we consider the Laplace transformation of this function. 
At first, we assume that p- m , p\ and P2 are independent variables, and at last we input 
Pm = Pi + P2- The Laplace transform ^, m) n(a;Pi,P2) is given by, 



a 



(4-163) 



where 



<t>f(a,U)= I ^trQU)e~ a - p = 2n^K 2 (am f )tr$U) : 



m. 



e -aq = 



cr 



a 



4>h{ot) = J 



d 3 k _ , 2-nm 

H K 1 (am H ). 



2E k a 

The above functions K v are the modified Bessel functions and a 
perform the group integration easily, 

21-2 f r 1 2 



(4-164) 

{tr(^)} ? 

(4-165) 
(4-166) 

= Then we can 



where 



dmrfaU^TfaU) {tr(^[/)} {tr(^C/)} -a 1 

= jtmQ I dUti^UM^U) {tr$U)} 2l ~ 2+21 (-a 2 )""* 
= {C (pi • a)(p 2 ■ a) + CM ■ p 2 )a 2 } a 2{l+m ~ 2 \ 

n — a W \m-i n X^v 1 * ~~ 



(4-167) 



C 1 = 2j2(-) m - t m C t - 



{2(Z + <-!)}! 



(4-168) 
(4-169) 



(l + i + m + i-l)l 

Finally, using the inverse Laplace transformation, we obtain the function Ii >m>n . We can 
evaluate ij jmjTl analytically in the case of the relativistic limit and the non-relativistic limit. 
In the relativistic limit, Ii >m . n is given by, 



T _ r ER 



ER 3i+2m+n-4 

,m,r 



(4-170) 



where 



c er = 2 r. 



(i) 



2l+m+n— 3 



m 
i=0 



m—i /~i 



{(3/ + 2m + n-3)!(3/ + 2m + n-5)!} 1 



(i+i+i)\(i+i-iy. 



{! + (/ + i)(3Z + 2m + n - 4)} .(4-171) 
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In the non-relativistic limit, we obtain 



where 



Il, m ,n = C^ n m 3 / 2 m 5 ™/ 2 m n J 2 (^ - M f + ^ m+n ^ w ■ s 1 / 4 , (4-172) 



Cf* n = | (el + \m + \ - lo) l} 1 • 2 il ^ m+n y 2 - 7 n 2l+3 ^/ 2 - 3 



in 



x EC")- M \^ ^ (l + i). (4-173) 



{2(i + «-l)} l 

In the general cases, we can perform the inverse Laplace transformation numerically, using 
the steepest descent approximation. We show the result in Fig. 19 for the case where n w = 75 
and rih = 0. In Fig. 19, the solid line denotes the numerical result by the steepest descent 
method. The dashed line and the dotted line denote the approximate result in the non- 
relativistic limit and extremely relativistic limit, respectively. We show cross sections for the 
various values of n w and nh in Fig. ^0|at ^=15.55 TeV. From this, we understand that the 
gauge boson plays an important role rather than the Higgs particle. In Fig. 21, we show the 
energy dependence of the cross section for n w =0, 20, 40, 75 and 100, and n h =0. We show 
the unitarity bound as a dashed line in Fig. 21. The approximation breaks in the region of 
the energy where the cross section overcomes the unitarity bound. In the region, the multi- 
instanton effect & (see section 2) or interaction between the gauge bosons in the final stated 
are important, which are not evaluated in our analysis. Since the asymptotic behavior of the 
valley instanton is different from the constrained instanton, a naive calculation yields that 
the interaction between the valley instanton and the anti valley instanton is different from 
that between the constrained instanton and the anti constrained instanton. Therefore, the 
analysis given in section 2 must be reconsidered. This is interesting problem, but remains 
as an open question for now. 



§5. Complex-time methods 

As we have seen so far in this paper, the imaginary-time path-integral formalism com- 
bined with the proper valley method has been successful for the treatment of the quantum 
tunneling phenomena. In spite of this, some basic as well as practical questions related to 
analytic continuation remain. How does the boundary condition matches that of the val- 
ley configurations? In estimating many-point Green functions, how does the external fields 
enter (affect) the imaginary time configuration? [This is crucial for the validity of the path- 
deformations in the valley calculation, as well as bounce calculations.] Some other points 
are elaborated by Boyanovsky, Willey and Holman. 0) 
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Fig. 19. The energy dependence of the cross section for n w = 75 and = 0. The solid line denotes 
the result by the saddle point approximation. The dashed line and the dotted line denote the 
approximate result in the non-relativistic limit and extremely relativistic limit, respectively. 
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Fig. 20. The cross section for the various n w and n/j at s 1//2 = 15.55 TeV. 



In quantum mechanics, the complex-time method has been studied by various authors. 
f|| , 0):il)S, IB >H)>il) It was argued that it allows semi-classical approximation for tunnel- 
ing phenomena and overcomes various problems associated with the pure-imaginary-time 
method. In it, one considers the analytic continuation of the time-integral in Fourier trans- 
form of the Feynman kernel. The existence of the tunneling leads to the existence of the 
complex saddle-points in the time-plane. It was claimed that semi-classical approximation 
is done by deforming the time-integral so that it goes through all such saddle-points. © 

This gives one hope that calculations in field theory may be improved by using the 
complex-time method. Indeed, Son and Rubakov0)-§)S adopted a complex-time method 
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n w =40 n w =15 n w =100 



Fig. 21. The energy dependence of the cross section for n w =0, 20, 40, 75 and 100, and n^=0. A 
dashed line denotes the unitarity bound. 

for the estimate of the cross section. They used periodic instanton-anti-instanton solution 
in imaginary time, and the initial Minkowski state overlaps with the instanton as a coherent 
state. 

Closer look, however, reveals many questions unanswered so far.© It is known that the 
complex-time plane is plagued with singularities and infinite number of saddle-points, which 
form lattice structure. Among saddle-points, there are two kinds, those that are solutions 
of the equation of motions (we shall call these "physical saddle-points" in this paper), and 
those that are not ("unphysical saddle-points"). How the path is deformed to avoid the 
singularities and to go through only the physical saddle-points is unknown. One simply 
assumes that all and only the physical saddle-points contribute. Even so, the weight of each 
saddle-point is a riddle. They are determined so that the result agrees with that of the WKB 
approximation. It is not known how it comes about from the path. 

In this section, we give the solution of these problems, based on a reduction formula 
in the number of the turning points of the path. In the first subsection, we give a brief 
overview of the saddle-point calculations in complex-time. There we elaborate some of the 
points briefly mentioned above. In order to solve the riddles pointed in this subsection, we 
analyze the ort honor mality of the WKB wavefunctions and construct the Green function 
from them in the subsection 5.2. This way, the arbitrariness in the path- integral formulation 
is avoided and the relevant Green functions are constructed from the first principle. Based 
on this construction of the Green functions, we derive the connection formula for the WKB 
wavefunction in the subsection 5.3, where it leads to the reduction formula for Green function. 
Expanding this formula, we obtain a series, which can be interpreted as a sum over the 
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physical saddle-points. The weights and the phases of saddle-points are determined from 
this formula. 



5.1. Saddle-point method 

We consider one-dimensional quantum mechanics with action; 



dt 



-x 2 - V(x) 



(5-1) 



The potential V(x) is assumed to be smooth enough to allow WKB approximation in asymp- 
totic regions I (x — > — oo) and II (x — > oo). The finite-time Green function (Feynman kernel) 
is defined by the following in the Heisenberg representation; 



G(x i ,xf,T) = {x { \e- iHT \x i ), 



(5-2) 



in terms of the Hamiltonian of the system, H. The advanced and retarded resolvents are, 
respectively; 

ro .._ 1 



G A (x u xf, E)=i j dTe i{E - iS)T G{x h x { ; T) = (x { 

J —oo 
roc 

G R {x h x f , E) = -i dTe^ E+ ^ T G(x u x { ; T) = {x u 

Jo hi + to — H 



(5-3) 



E-i5-H 1 17 ' 

Xi) , (5-4) 



where real-positive 5 is introduced to guarantee the convergence of the integrals. The poles of 
these Green functions come from the bound states, and the cut from the continuous spectra. 
Let us first look at the saddle-point approximation for the path integral, 

rx(T)=x f 
/x(0)=Xi 

One obtains the following from the saddle-point method; 



rx(T)=xt 

G(x u x { ;T) = / Vxe lS . 

Jx(0)=X; 



(5-5) 



G(x h xf,T) = J2 



X C 1 



2tt x cl (0)x cl (T) 
i -d 2 S cl /dT 2 



'1/2 



(5-6) 



where x c \ is the solution of the classical equation of motion and S c \ is the action evaluated 
at the classical solution. This approximation is valid as long as both x\ and Xf are far from 
the turning points where classical velocities vanish. Substituting (5-6) into the retarded 
resolvent in (5-3), (5-4), one arrives at the following expression; 



G R {x u x f -E) 



roc 

: / dry 



2tt x cl (0)x cl (T) 
-d 2 S cl /dT 2 



'1/2 



J(ET+S cl ) 



(5-7) 



In the above, we absorbed the infinitesimal imaginary convergence factor (iS) in E. 
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In calculating (5-7), one hopes to apply the saddle-point approximation to the T-integral. 
The saddle-point condition (the stable-phase condition, to be more exact) is then, 

dS cl 1 



E 



Ki + v(x cl ) 



(5- 



dT 2 

In other words, the solution of the equation of motion x c \(t) has to have the energy E 



Let us consider a situation depicted in Fig. ^2[ Both the initial point X\ and the final 
point Xf lie in the forbidden region. There is no real-time solution of the equation of motion. 
Thus no semi-classical approximation is possible for the real-time formalism. 

If one assumes that the expression (5-6) gives correct analytic continuation of the Green 
function in the complex T-plane when complex-T solutions of the equation of motion exist, 
one could go to the complex T-plane and apply the saddle-point method. One then arrives 
at the following expression; 



G(xi,x { ',E) 



1^1(0)^1(^)1 



-1/2 



w(T s ) p(T.) e m ^\ 



(5-9) 



where T s is given by (5-8). Here W C \(T B ) is the WKB phase, w(T s ) is weight and p(T s ) 
is phase for the corresponding saddle-point. The original contour is deformed to pass the 
saddle-points by a series of steepest descent paths from the origin T = to the positive 
infinity on the real axes (Re T = 00, Im T = 0). The weight of a saddle-point, w(x c \), is 
determined by how the contour crosses the saddle-point: If the contour does not cross the 
saddle-point, the corresponding weight is zero. If the contour crosses it along the steepest 
descent direction, the weight is 1. In case the contour reaches a sub-dominant saddle-point 
from the direction orthogonal to the steepest descent and leaves it along the steepest descent 
direction, the weight is 1/2 (as in the calculation of the false vacuum decay). The phase of 
a saddle-point, p(T s ) comes from the square root in (5-7). 

Carlitz and Nicole applied this method to linear potential, quadratic well and quadratic 
barrier, all of which are exactly solvable. Using the exact expression for the Green function 



Von 




Fig. 22. A potential V(X) with asymptotic regions I and II. 
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(5-7), they found that the saddle-points explained above, namely the ones associated with 
the complex-time solution, indeed lead to saddle-points in the T-plain [we shall call these 
"physical saddle-points"]. They, however found more: there are saddle-points that are not 
associated with any solutions of equations of motion (unphysical ones), and also singularities 
due to the periodicities. Deforming the T-integration path to avoid singularities, they found 
that the path indeed goes through only the physical saddle-points. As a result, the weights 
and phases in (5-9) are determined, which leads to the correct WKB result. 

In the case of a double-well potential, the deformation of the integration contour is not 
specified, which is understandable in view of the fact that the analytic structure of (5-2) in 
complex T-plain is fairly complicated. The claim made in the literatures is simply that the 
contour passes all saddle-points that correspond to classical trajectories with weights that 
are extracted from the above calculation. 

The problem lies in the fact that we do not know the contour which passes all physical 
saddle-points that are distributed in complex-time plane. Because of this, we cannot deter- 
mine the weights of the saddle-points. In the case when x- x and xt are in the forbidden region 



on the right side of the double- well potential as depicted in Fig. g2], saddle-points in complex 
T-plane are given by 

T s = ±T(b 2 , xi) ± T(b 2 , x { ) + lT(a u h) + mT{b u a 2 ) + nT(a 2 , b 2 ) (5-10) 
l,m,n = 0, ±1, ±2, ... , (5-11) 

where 

ry cIt' 

T(x,y) = 2 (5-12) 
J' p{E-V{x')) 

Notice that T(x, y) is pure imaginary if (x, y) is in a forbidden region. These saddle-points 
T s are shown as solid circles in Fig. The unphysical saddle-points are also shown as open 
circles in the same figure. We know of no contour which passes only the physical ones and 
avoids the unphysical ones. Thus so far the complex-time method has not been worked out 
for the double-well and more complicated potentials. 

5.2. Green function and the WKB wavefunctions 

We shall examine the Green functions by using the complete orthonormal set of the 
eigenfunctions {l^n)} of H. For example, the retarded Green function in (5-3), (5-4) is 



written as 



G ix u x u E)-^ E + i5 _ Xn , (5 13) 
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Fig. 23. Positions of saddle-points in complex T-plane. 

where the sum over the complete set is actually made of the integration over the continuum 
spectrum and sum over the discrete bound states. 

In order to study the Green functions in the asymptotic region, we use the WKB ap- 
proximation for the wavefunctions ip(x) = (x\ip). For a state with the eigenvalue A in the 
continuum spectrum, the second-order WKB approximation yields the following; 

(Ax^C p(x ' )dx ' + Bxe-'f* 1 for x el, 



2irp(x) 



C x e 



+ D x e 



i f X p(x')dx' 



for xE II, 



where p(x) = y 2(A — V(x)). We choose the limits x\ and X2 of the integrations to be the 
nearest turning points for definiteness. Among the coefficients A, B, C, and D, only two 
are independent. Their inter-relation is linear due to the superposition principle and can be 
written as follows 0; 



(5-15) 



Ax ] = S(X) I Cx 

The 2x2 matrix S(X) is determined by the shape of the potential V(x) in the intermediate 
region between x\ and x^- The flux conservation law is satisfied for the Schrodinger equation; 



- • - 

dx^ ' 



2 \ dx dx 

Substituting the WKB expression (5-14) into the above, we find that 

3 = ^{\A\ 2 -\B\ 2 ) = ^{\Cf-W). 
This results in the following relation between the matrix elements of 5(A); 



\S\i\ — \Szi\ — 1, \S\2\ 

S\iS\2 — 5*2! 5*22 = 0. 



\S- 



22 



(5-16) 

(5-17) 

(5-18) 
(5-19) 
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As a result, the matrix 5(A) is parameterized by three real functions a(\), f3(\), and p(A) 
as follows; 

/ e M cosh p e l(i sinh p \ ^ 
I e~ t/3 sinh p e _M * cosh p y 

For a given eigenvalue A, there are two independent eigenfunctions. In order to con- 
struct them, we look at the inner-product of the eigenfunction ip\(x) with coefficients 
(A\, B\,C\, D\) and another eigenfunction ip\>(x) with (Ay, By, Cy, D\>). Integration in 
the asymptotic regions determines the coefficient of the delta function 5(X — A') completely, 
as the delta function can come only from the infinite integrations. [In order to calculate the 
finite terms we need to solve the Schrodinger equation completely. But this is not necessary 
for the current purpose.] In fact, integration in the region II yields, 



f°° ~ 1 f°° dx 

/ dx *l>\(xm>(x) = — / , 

Jx 2 l-K Jx 2 ^p(x)p'{ 



X 



'(x) 

n*n jf*Jp( x ')-p'( x ')) dx ' , n* n - i f*(p( x ')-p'(. x ')) dx ' 



V 



+ C* x Dye I > ix ' )+P ' {x ' ))dx ' + Dl Cye- l ^ ip{x,)+p,{x ' ))dx ' I , (5-21) 



As we are interested in the singularity of the above for A = A', we expand the exponent using 
the following for A ~ A', 

P(x')- P '(x')c { -^1. (5-22) 

By using a new coordinate y defined by dy = dx/p(x), we find the contribution from the 
first two terms to the delta function; 

rco _ 1 roc , „ . 

/ dx iPl(x)$ x .{x) = — dy (C* x C x e t{x - x ' ){y - y2) + D* x D x e-^ x - x '^- y2 A + ... 

Jx 2 27T Jy 2 V ' 

= \(ClC x + DlD x )6(\-\') + ..., (5-23) 

where we have neglected all the finite terms. Doing the similar calculation for the region I, 
we find 

dxrx(x)^y(x) = - (A* X A X + B* X B X + C* X C X + D* X D X ) 5(X - A'). (5-24) 
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Using the above result, we choose our orthonormal eigenfunctions for a given A as the 
following; 



(A X ,B X ,C X ,D X )^ 



, 0, 1, -e i(a-/?) taohp) , 



v cosh p 

(A X ,B X ,C X ,D X )V = (e^tanhp, 1, 0, 



cosh p 



(5-25) 
(5-26) 



[The above corresponds to two incoming states. There are of course other choices, such as 
stationary states, but all of those reads to the same result for the Feynman kernel.] 

Let us first look at the case when both x\ and Xf are in an allowed region in II. From 
(5-25) and (5-26), we find that 



(i)* , 



i=l,2 



X 



2n^p(xi)p(x f ) 

p P (x')dx> i(r P (x')dx'+nP(x')dx') Ha-fj). , , , \ 

Jx i - e V*2 J *2 ) e na P) tanhp + (c.c.) 



(5-27) 

(5-28) 
(5-29) 



The complex conjugate part can be understood as the same function below the cut on the 
real-axis of the complex A-plane (see Fig. 24). This is guaranteed by the fact that as A 
moves below the cut, the role of the coefficients A and B is exchanged and likewise for C 
and D, and as a result, the phases a and (3 change their signs, while p does not. Therefore 
the sum over the continuous spectra in the Feynman kernel (5-13) is written as the following 
A-integration; 



+ iS - A 



p(xi)p(x{) 



E + iS 

X 



C 2 Ci 



Fig. 24. The complex-plane of the eigenvalue A. The cut on the real-axis is the continuous spectra 
and the poles are the bound states. 
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x < 



i f Xi p(x')dx' 

g Jx 2 

i P f p(x')dx' 

g Jx 2 



-i f X{ p(x')dx' i f Xf p(x')dx' ;( n _a\ , , 

e Jx 2 — e Jx 2 e [ ^tanhp 

—i f Xl p(x')dx' i f Xl p(x')dx' j( n —ff\ , i 

e Jx 2 _ g Jx 2 g ^ ^tanho 



for Xi > Xf, 
for Xf > Xj. 



(5-30) 



The integration contour C in the complex A-plane is as in Fig. 24. In (5-30), we have chosen 
the integrand so that it converges for |A| — > oo, which we need later. As for the discrete 
spectrum, it is determined by the absence of the diverging behavior of the wavefunction, 
which translates into the condition that for C = 0, A = 0. From (5-20), we thus find that 
sinh p = when A is equal to a discreet eigenvalue of H. This allows us to write the sum 
over the discreet spectra as pole integrations of the second term in (5-30) depicted in Fig. 24. 
Connecting all the contours and closing it at the infinity, we find that the whole contour 
enclose the pole at A = E + i5. Thus we end up with the expression 

i 



G R (x h xf,E) = 



p(Xi)p(X{) 

f Xi p(x')dx' ( -i t Xf p(x')dx' .„ i f Xf p(x')dx'\ £ 

J *2 FK (e J *2 y ' +tRe Jx 2 for x; > x f , 



x 



i f Xl p(x')dx' ( —i j Xf p(x')dx' . „ i f Xf p(x')dx'y . 

e Jx 2 ' (e Jx 2 py ' + t R e Jx 2 ^ > y ) forx f >Xi, 



(5-31) 
(5-32) 



where we defined the "reflection coefficient" R by the following, 

R 



i^i =ie i( - a -^ tanhp. 



(5-33) 



'22 



We note that this expression (5-32) is valid for a general value of E with proper analytic 
continuation of the coefficients. Therefore, if any end-points are in a forbidden region, simple 
analytic continuation of (5-32) is appropriate. We shall use this result in the next subsection. 

5.3. Reduction formula and its expansion 

In this subsection, we derive a formula for evaluating the Green function in a system 
with an arbitrary potential to which we can apply the WKB approximation. We consider a 
general potential which has arbitrary number of wells, and show that the Green function is 
given by summing up all contributions of classical paths. 

Let us consider the case when x- x and Xf lie in the forbidden region on the same sides of the 
wells. The potential which has n wells is depicted in Fig. 25. By the analytical continuation 
of (5-32) in E, we find the retarded resolvent in the following form; 

G R (x u x { ;E) = -\p( Xi )p( Xi )\-V> e - Ai (e A f +iR n e~ A f) (5-34) 

where 



f x i,f 

A j = / dx\p(x)\. 



(5-35) 
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In order to specify the fact that this expression is for n wells, we attach the subscript (n) to 
the coefficients hereafter. 

Due to the existence of the intermediate region (fe n _i, a„), the matrix S^ n \ which connects 
the regions I and II, can be written in terms of the matrices S**" -1 ) that connects regions I 
and (& n _i, a n ) and S for (6 n _i, a n ) and II. If we apply the linear WKB connection formula for 
the latter region, which we can obtain from the saddle-point method for a linear potential, 
we obtain the following; 

-An-l 



x 



where W n and A 



= S {n - l \E) 
n _i are defined by, 




A»-i 



b n -i 



dx\p(x)\, 
dx\p(x)\. 



Therefore, we find the relations; 



S- 



1 



(n) 
21 

o( n ) _ 

J 22 ~ J ' 



s. 



(n-l)-A v 



21 



sin W„ 



(n-l) 
21 



e^"" 1 sin W n + 2S%~ 1] e^"" 1 cos W n . 



From (5-33), (5-41) and (5-42), we find that R n can be written as the following; 



Rn — 



sin W„ 



-2A r , 



l R n -i cosW n 



2 cos W n - \e- 2A ^R n _ x sin W n 



(5-36) 
(5-37) 
(5-38) 

(5-39) 
(5-40) 

(5-41) 
(5-42) 

(5-43) 



Voc) 




E 



b l a 2 b 2 4,-1 b n . x a„ b n x f x t 



Fig. 25. A potential with n wells we consider. The turning points are denoted by a« and 6j, for 
left and right of the i-th well, respectively (i = 1 ~ n). 
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We find it most convenient to rewrite this to the following two expressions; 

1 1 — i? e 2iWn 

1 y (5 . 44) 

2 l + R n e 2iW ^ V 7 
1 - \R n -ie~ 2A ^ 

Defined in this manner, the function R n -\ corresponds to the reflection amplitude at the 
turning point b n -i- 

In order to see the correspondence between these expressions and the saddle-point 
method, let us expand (5-44) in an infinite series as follows; 

iR n = % - + (-iR n ) e 2m - + (-iR n ) 2 (-i) + ... . (5-46) 

The convergence of this series is guaranteed by the implicit factor iS. The retarded resolvent 
is then, 

G R (x i ,xf,E) = -\p(x i )p(x t )\- 1 ' 2 



x 



+ |1 + (_^ n ) e 2m ~ + (-iR r ) 2 H) e w " + ...} e~^ +A ^ 



(5-47) 



This is the expression to be compared to the saddle-point expression (5-7). If we choose T 
to be negative imaginary, —ir, the factor in the exponent is, 

i(ET + S cl ) = Et - S E = - J dr (^j =-Jdx\p(x)\. (5-48) 

Therefore, the two factors A it f is equal to the above quantity for the path from x it f to 
the turning point b n . The first term in (5-47) corresponds the contribution of the pure- 
imaginary-time path that starts from x\ and reaches Xf directly. The factor e~( Ai+Af ^ in 
the rest of the terms is for the path from x\ to the turning point b n and then from b n to 
Xf. The expansion of iR n is understood as contributions of the paths that oscillate in the 
allowed region (a n ,b n ). Various factors have unique interpretations as factors coming from 
the turning points. This is most conveniently depicted in Fig. 26. Similarly to (5-46), the 
expression (5-45) is expanded as the following, 

-iRn = -t + (iRn-!) e~ 2A -i + (tRn-i) 2 ^) e^ A ^ + ... . (5-49) 
The corresponding diagrams are illustrated in Fig. 27. 
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Fig. 26. Diagrammatic representation of the expansion of iR n . The information of the region left 
of a n is contained in —iR n . 




dn-\ bn-\ a n bn X f X t 




Fig. 27. Diagrammatic representation of the expansion of —iR n . The each path corresponds to a 
group of paths with different number of oscillation in the forbidden region (b n _i, a n ). 

We have so far derived the expressions that results from the WKB approximation for the 
n-th well. This procedure can be applied recursively for the rest of the wells. At the end, 
we are left with Ri, which is given by 

iR 1 = \ tan W 1 = l - + {-%) e~ Wl + H) 3 e" Wl + ... , (5-50) 

where W\ is denned by (5-39) by replacing n with 1. Therefore we find that R\ is expressed 
as a sum of the contribution of the classical paths which evolve in the first well. 

Combining (5-47), (5-49), and (5-50), we find all complex-time paths are included in the 
expression of the resolvent G R . Therefore, the resulting Green function in a n-well potential 
is given by 

G R (x uXt ;E) = - IpOrO^f )r V2 E /We WMiJ!) . (5-51) 
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where f(x c \) is determined by the number of reflections and transmissions which the classical 
path contains, and W{x\, xf, E) is determined by what wells the path crosses and what 
barriers it tunnels through. The rules for calculating f(x c \) is the following: When a path 
has a reflection in the allowed region, it obtains —%. In the case of a reflection in the forbidden 
region, it gets i/2. 

Let us next examine the case when x;, X{ are on opposite sides of the wells. Applying 
the analysis similar to the one in the previous section , we find that the Green function is 
expressed as follows; 

G R {x hXi -E) = -WxMxtT^ne-*-*' ■ (5-52) 
where T n is the (analytically-continued) transmission amplitude 

T n = * . (5-53) 

^22 

Just as the previous case, we can express T n in terms of T n _i as in the following; 

T n = — t — tn ( 5 ' 54 ) 

V^-* sin W n + 2S&-%a»-i cos W n 

rp -A n -!+iW„ 

(5-55) 



1 + \R n ^e- 2A ^ + e 2iW - (l - \R n ^e- 2A ^) 

= i ^ T n _ v (5-56) 

1 + |i? n _ ie - 2 ^-i 1 + R n e 2 ^ n 1 ; 

In the above, we have used (5-33), (5-42) and (5-45). Let us show that T n consists of all 
contributions of classical paths. The second factor in (5-56) is expanded in the following 

way, 

— 5 2lW =e^ + H)(-^ re )e 3 ^ + H) 2 H^) 2 e 5 ^ + ... , (5-57) 
1 + R n e ztWn 

which corresponds to the diagrams in Fig. 28 (a). The first factor in (5-56) is expanded as, 

^—A n -i ■ ■ 

t n / . „ n Q A . ,6 



= e^- 1 + 0(iRn-i) e- 3A "-> + Q 2 (*i?„-i) 2 e- 5A "-> + .(5,5? 



1 + \R n ^ e-^-i " ' V^"""- 1 ^ ' v 2< 

which is illustrated in Fig. 28 (b). As the previous case, repeating this procedure, we reach 
Ti, which is given by 

Tx = — — j— — (5-59) 
2 cos W^i 

= e^ 1 + {-if e 3iWl + H) 4 e 5lWl + ... . (5-60) 
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(a) 




Xf Cln-\ bi-l Cln bn X- t 

Fig. 28. Classical paths traversing in the n-th well and tunneling through the (n — l)-th barrier. 




Therefore we conclude that T n consists of all contributions of classical paths. 

The above analysis can be applied to other cases with various locations of x- x and Xf. 
Thus we verify that thus constructed resolvent can be always interpreted as sum over the 
complex-time classical paths. 

This formalism is valid for a scattering process and metastable system (Fig. |29[) as well as 
for a stationary system. It is apparent from this derivation that this complex-time method 
reproduces the result of the WKB approximation with the linear connection formula. The 
potential we consider is given in Fig. (b). The turning points are labeled in the figure 
by a 1; bi, and a 2 . First we consider the case when X\ is in the well and Xf is outside of the 
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barrier. The analysis similar to the above yields the Green function as in the following; 

f e -iWi _ ie iWi\ e iWi-A!+iW f 

G R (xi,xf,E) = -iMxMxAr 112 ^ = '- = , (5-61) 

where R is given by 

1 - -e~ 2Al 

= -i + V 2 ^ + (^) 3 e- 4 ^ + ... . (5-63) 

From (5-61) and (5-63), we find again that the resolvent is equal to the sum over the physical 
saddle-points. We find that the poles of the Green function are determined by the following; 

1 + e 2 ^ + l -e~ 2A ' (l - e Wl ) = 0. (5-64) 

Let us solve this equation iteratively. Then we obtain 

Wi{E) = ( n + ^) 7r ~ \ e ~ 2Al - ( 5 ' 65 ) 
From (5-65), we find an imaginary part of the energy eigenvalue 

lmE n = % ——e' 2A \ (5-66) 

2T(E n ) ' v ; 

where T(E n ) is the period of the classical path between the turning points a±, hi 

T{E) = 2 (5-67) 
J <* y/2 (E - V{x)) 



We find the decay rate of this metastable system. We note that the factor 1/2 in (5-66) 
comes from the weight 1/2™ of the saddle-points with n reflections in the forbidden region. 



Let us apply this method to a scattering process against a potential barrier in Fig. g9 
(a). When x\ and xt are separated by the barrier, the Green function is given by 

G R (x i} x f ,E) = -i|p(z i )p(a; f )|- 1/a fe i ( Wi+w » 

= -i\p(xMxf)\~ 1/2 Y, ei{Wi+WfM2k ~ 1)Al (-) ( 5 - 68 ) 

k=l 

where we have used the WKB expression for the transmission coefficient T; 

f = ttt^- (5 ' 69) 
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When x\ and X{ are on the same side, the Green function is given by 



G R (x h x t ; E) = -i\p( Xi )p(xf)\- 1/2 j Wi [e~ iW f - iRe iW f' 

■ oo 

= -i\p( Xi )p(xf)r 1/2 e m e- lW f + % -e lW f + £ e lW f- 2kA ^-t) 

^ i i 



fc=i 



2k- 1 



(5-70) 



where we have used the expression for the reflection amplitude, 




(5-71) 



The coefficients T and R satisfy the unitarity condition 



T| 2 + \R 



1. 



(5-72) 



We again confirm the validity of the sum over the physical saddle-points. The weight l/2 n 
is crucial for (5-72). 

We have given the re-formulation of the complex-time method. In quantum mechanics, 
using the connection formula for the wavef unctions, we constructed the reduction formula 
in the number of the turning points for the Green function. This yields series expansions, 
which can be understood as sum over in the classical complex-time trajectories. This is 
understood as a sum over the physical saddle-points with specific weights and phases in the 
path-integral method. This shows the validity of the method proposed before, in the context 
of the path-integral method. We confirmed that this method yields results identical with 
that of the WKB approximation. Thus our construction gives solid basis for the starting 
point of the complex-time method. 

We must elaborate this complex-time formalism so that it could be applied to field theory 
and two-particle scattering process. Thus somewhat different analysis may be required. One 
of them could be the combination of the valley method and the current complex-time method. 
The valley method can be used to identify the imaginary-time tunneling paths, which are 
converted to collective coordinates. Thus, while the incoming particle are expressed in real- 
time expressions, the tunneling part may be obtained as the imaginary part of the complex- 
time development along the valley trajectory. There are rather interesting possibilities along 
this line, which should be further investigated. 



The tunneling phenomenon has some subtle relation to the large order behavior of the 
perturbative series of the theory. Dyson argued© that the perturbative series in QED is 



6. Asymptotic analysis of the perturbation theory 
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an asymptotic series. Later, more accurate estimates of the large order behavior of the 
perturbative series were carried out in many models.© Although asymptotic series are di- 
vergent, it allows summation techniques, such as Borel summation and Pade approximation. 
In theories with tunneling phenomena, the perturbative series, however, are generally not 
Borel-summable. In such a case, the Borel integral encounters singularities, which are related 
to the existence of the nonperturbative tunneling phenomena. Therefore, only after the per- 
turbative and nonperturbative effects are combined the physical quantities are well-defined. 

The consideration of the valley, described in detail in this paper, provides an insight into 
this interplay between the perturbative and non-perturbative contributions to path-integral: 
In the context of the the high energy behavior of the B and L violating amplitude considered 
in section 2, we have noted that the contribution of the configuration which has an instanton 
and an anti-instanton in small separation compared with their size become important. Since 
this configuration merges into the vacuum continuously, it implies that the perturbative and 
non-perturbative contributions transmute into each other in this configuration. This situa- 
tion casts a doubt against our intuitive understanding that there must be a clear separation 
between the two. For the full understanding of the tunneling phenomena, we consider it 
necessary to understand this situation. Approaches in this section provide some possible 
direction for progress. In the first subsection, we show that the valley calculation provides 
an easy way to evaluate the large order behavior of the perturbative expansion in a model in 
quantum mechanics. 0* In the second subsection, the asympton theory, which is an analysis 
of the perturbative functional, not just the coefficients, is presented. 

6.1. Borel function in the valley method 

We consider the symmetric double-well potential 



and evaluate the transition amplitude (0|e HT \0) by the path- integral Eq. (|4-2|) . We are 
interested in the perturbative expansion of 



(H = (1/2) [p 2 + 4> 2 } the free part of the Hamiltonian) or the energy shift E(g 2 ) of the ground 
state due to the interaction, which is evaluated by 



The time T is taken sufficiently large, compared with the inverse of the excitation energies. 






(6-2) 




(6-3) 
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The perturbative series of Z(g 2 ) in the present model is a power series in g* 



Z peit (g 2 )=l + J2Zn9 2n . 



(64) 



n=l 



The coefficient Z n are evaluated by expanding the action in powers of g, 

e- s ^ = Y.9 n F n [^ S[hg] = (* dr\U$f + V(hg) 

Jo 12 \ ' 

and path-integrating each term, 



JV(fiF< 



2n\ 



in the background 0b = 0, where So is the action for the harmonic oscillator 

/■T 



So 



dr- 



(</>) + <P' 



(6-5) 



(6-6) 



(6-7) 



Note that the functional F n [(p] with odd n has odd powers of <fi and gives no contribu- 
tion. Although Z n are well-defined objects, the infinite series does not converge; 
§,§),§),§),§),§) it ig known tQ be asymptotic to Z(g 2 ). 

The Borel summation method is a powerful one to analyze the asymptotic series Eq. ( |64|) . 
To apply the method, we first calculate the Borel function B(t) defined by a series 



Z„ 



(6- 



Z(g 2 ) is recovered, at least formally, by the integral 



Z(g 2 )= dte- t B(tg 2 ) = ^ dte~ x ^ ? Bit) 
Jo q 2 Jo 



since 



dt e~H 



t±n 



(6-9) 



(640) 



Note the series (|6-8|) has wider region of convergence in the complex t-plane than the series 
(|6-4| ) has in g 2 -plane. This may give a chance to get a meaningful answer even from the 
non-convergent series. 

We intend to calculate the Borel function of the present model not by the series expansion 
( |6-8| ), but in a rather direct way using the valley method. It is well-known that the present 
model possesses the instanton (I) and anti-instanton (I), the solutions of the Euclidean 
equation of motion that travel between = and (f) = 1/g. Applying the proper valley 
equation, we have further found there exist the valley trajectory that starts as a small 
fluctuation at the vacuum and reaches to the 14 pair at the end. 0) We first approximate 
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Z(g 2 ) to an integral over the valley trajectory 0y(a), where a can be any variable that 
parameterize the valley In another word, we adopt a line 4>v{oe) in the functional space as 
the background for calculating Z(g 2 ), instead of using a single point 0b as in the perturbative 
calculation. The explicit expression as the collective coordinate integral is obtained by 
the Faddeev-Popov (FP) method as explained in section 3. To avoid the cumbersome g 
dependence in the background valley </y(a), we rescale S and <fi as 



S[<f>;g]^g 2 S[<f>). 



(6-11) 
(6-12) 



With this transformation, g drops out from the action and 0y(a) becomes g invariant. We 
now obtain 



dae-sWa)]/? 



dS[</>v((x)] 



da 



8S 



Z. 



(6-13) 



where 



Z 



Af 



27rc/)(0|e-^ T |0) 

g(ip ■ D(a) ■ d4>y(ot) / da 



J [c^j'exp (-^(<f- D(a) ■ <p) - S?Jp; 9) 



1 - 



dS[(py(a)]/da 



(6-14) 



In Eq. (6-14), <p = g(<fi — 0y(a)) represents the fluctuations around <pv, [dtp]' is the functional 
measure except the mode parallel to 8SJ 8<j) that was removed by the FP method, 1 / V2ng 
comes from this elimination, D(a) represents 5 2 S/5<p v (a) and the abbreviation like <p-D(a)-ip 
stands for 

i-T A 2 .Q 

— <p(T)ip(l J ), 



drdr' — — — 

o 5(f)v{a)(T)5(J)v{a)(T' 



S' mt is the remaining action after the value at <j)y and the quadratic term (pDtp are subtracted. 

The trick of extracting the Borel function from Eq. (|6-13|) is to choose the parameter a 
as the value t of the action of the valley, so that t = S[0y(£)]. Then Eq. (|6-13| ) is written as 



[ Scl dte-^ 2 


5S 


8<f>v(t) 





z. 



(6-15) 



where S c i — 1/3 is the twice of the instanton action which the valley action cannot exceed. 
Note the resemblance of this expression to Eq. The integrand 



B v it) 



8S 



-i 



Z 



8Mt) 

is the approximation for the Borel function in terms of the valley method. 



(6-16) 
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Although By is an approximation, it suffices to predict the large order behavior of the 
perturbative expansion. The bottom line is that the large order behavior is controlled by 
singularities of the Borel function in the complex t-plane and that the valley approximation 
By also reproduces them. The valley approaches asymptotically to the well-separated I-I 
pair <fiii(R) with R their separation. The action t behaves 

t ~ S d - const x e~ R (6-17) 

at There exists the normalized quasi-zero mode 0r that causes variation of R, i.e. 
0r = (l/r])(d(j)ii/dR), where 77 is the Jacobian between R and 0r. The functional gradient 
is proportional to the quasi-zero mode according to the definition of our valley equation, 
0r = (5S/54>y)/\5S/5(f)y\. Since (5S/5(py) ■ (d^i/dR) = dt/dR, we can conclude 



5S 



SMt) 



12(S d -t) 



(6-18) 



where we have used r] = yl/12. The calculation of Z is also simple at the well-separated 
I-I pair; to the leading order in g 2 it is approximated to the product of the square of the 
single instanton determinant and the Jacobian for the overall translation, 

1 12 

Z = -—^=T. (6-19) 

2th? 2 v/3 1 ; 

(The calculations of the determinant and the Jacobian are found in Ref. §).) Combining 
these, we obtain 

B v^) - —2^r—+- ( 6 - 2 °) 

ng 2 Sci - t 

This simple pole controls the large order behavior of the perturbative expansion of Z(g 2 ). 
In fact, inserting its expansion ~ (T/7ig 2 )J2n{Sci)~ n ^ 1 t n into ( |6-15| ) and integrating it over 
to infinity (neglecting the upper bound), we obtain 

Z n ~ -S A - n - l n\, (6-21) 

7T 

or, for coefficients of perturbative expansion of E(g 2 ) in powers of g 2 , 

E n ~ --S^nl (6-22) 

7T 

which coincides with the known results.© 

The singularity of By(t) at t = S c \ indicates the non Borel-summability of the expansion. 
There will be an ambiguity depending on how we turn aside the pole when we perform a 



contour integral corresponding to Eq. ( |6-9| ) if we only know B(t) from the perturbative ex- 
pansion. The advantage of using the valley method is that we now know the ambiguity is, 
at least in the present model, an artifact: the integral has the upper limit at t — S c \- This 
understanding may provide us to construct a specific prescription to get unambiguous predic- 
tions from the perturbative expansion by properly taking into account the non-perturbative 
contributions in a class of models which is known to be non Borel-summable.0 

6.2. Asymptons 

In this subsection, we look into the behavior of the perturbative functional , F n [<f>] in 
Eq. (|6-5[), in the functional space of path- integral and discuss it is dominated by a series 
of configurations, which we name "asymptons". We will show their contribution indeed 
reproduces the large order behavior of the perturbative series, which shows the fictitious 
nature of the non Borel-summability and provides us with the insight into the interplay 
between the perturbative and the non-perturbative effects. 

In the first part of this subsection we present a general formalism for the analysis of the 
asymptons. Then a simple double-well quantum mechanical model are treated in the next 
part. The SU(2) gauge field theory is studied in the third part. Then the analysis of the 
asympton on the valley line is carried out. This subsection ends with the discussion on the 
consequences of the asymptons and on the summing methods. 

6.2.1. General formalism for the asympton analysis 

Let us use z n for each coefficient in the perturbative expansion of a given model (not 
necessarily the symmetric double- well as in the previous section), 

_ fV^FM 

We note that the coefficients z n are known to grow as a factorial of n for large n in a wide 
range of models while the expansion (|6-5f) is convergent for a given configuration as long as 
the integrals J dr(0) m are convergent. Our aim here is to understand how this can happen. 
For this purpose we study the behavior of the functional F n [<j)} itself. The asympton theory 
answers this question. 

The expression for F n [(j)] is complicated in general. In the limit n — > oo, however, a 
compact expression can be obtained. The functional F n [0] is expressed as a contour integral 
in the complex (/-plane as follows; 



See Ref. |47j) for the detailed discussion on this point. 
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where 



S[<p,g] = S[(j),g) +nhgg. 



(6-25) 



We apply the saddle-point approximation to the g-integral ( 6-24 ) valid for large n. 

In this subsection we deal with theories whose action is a second order polynomial of g, 
which we express as S[4>, g] = cq[</>] — gci[4>] + g 2 c^)\. The saddle points of S[<f>, g] are given 
by the following; 

8nc2 



g± 



Cl 
4c 2 



[1±VD), D = l 



r 2 
c l 



(6-26) 



In order to examine the behavior of the perturbative functional, we divide the functional 
space of 4> into two regions; the one with D > and the one of D < 0, for the the expressions 
of F n [<f>] obtained in the saddle-point approximation are different depending on the sign of D. 
As was found in Ref. [51]), [52]), the region of negative D is of importance. We will concentrate 
on this case hereafter. 

When D < the saddle points are complex, which we denote by g± as follows; 



9± = \g\ e 



±i0 



77 



\g\ 



2V C ° S ^ 



Cl 



y/8nc 2 

where 9 G [0, 7r/2]. We obtain the following expressions at the saddle points g±; 
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c -n 
2n 



,±2id 



\g\ 



sin 9 (sin 9 ± icosi 



(6-27) 

(6-28) 
(6-29) 



Since Re(S^) > at both saddle points, we choose the integration contour to go through both 
saddle point along the real axes. The contour integration yields the following expression; 



-Am 



exp(-S[0,cj_]) 



S'Lgl 



(6-30) 



We stress here that g_ is a functional of </> and a function of n, defined by (6-27). 

Let us look for peaks of F n [(f>] in order to understand its the large order behavior. There- 
fore we must solve the following equation: 

5F n [0] dF n [<p] 5g- , 5F n 



6<t> 



dg- 



0. 



(6-31) 



At the leading order of H (which we have neglected to write), the g_ derivative in the first 
term acts on S and is zero due to the saddle-point condition. At the same order, we get the 
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following expression by using (6-30); 



Im 



SS[<p,g. 



(6-32) 



where 



a = are 



S[<t>,g- 



S-9 2 - 



9- 



1 • n ^ 9 * 
- - sin 20) + - + -. 

2 ; 2 4 



(6-33) 



This expression (6-32) is similar to the ordinary equation of motion except for the fact that 
the "coupling constant" g_ and the phase o are functionals of <fi. The solutions can be 
identified by first solving the equation of motion (6-32) for arbitrary constant g>_ and a and 
then solving the self-consistent equation (6-27) for \g\ and 9. 

6.2.2. Asympton for the double-well quantum mechanics 

We first consider the quantum mechanics of the double-well potential, whose action and 
perturbative functionals are given by Eqs. fl6-5|) and (6.1). Then the equation (6-32), which 
determine the peaks of -F n [0], can be written as follows by using the scaled variable ip = \g\<f>] 



d 2 (p dVes(<p) 



+ 



0. 



dr 2 dip 

The "effective potential" V e s((p) in the above is defined by the following; 

Vcs(p) = -<p 2 - K 3 p 3 + i/t 4 y? 4 , 



(6-34) 



where 



«3 



K4 



sm[o~ — a ) 

sin a 
sin(o- - 29) 

sincr 



(6-35) 

(6-36) 
(6-37) 



play the role of the "coupling constants" . The self-consistent equation (6-27) is now rewritten 
by using the scaled variable as follows; 



' V n ' 
R 3 



cos 9 



2R, 



(6-38) 
(6-39) 
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where 



R m = J drcp m . (6-40) 

We note that the effective equation of motion (6-34) is independent of \g\; it has only two 
free parameters, n and 9 (or alternatively n 3 and k 4 ). Consequently, the only self-consistent 
equation we need to solve is Eq. (6-39), which determines the value of 9 for a given n; 
Eq. (6-38) simply determines the value of \g\, which will be used for the evaluation of F n [0]. 

The effective equation (6-34) possess a bounce-like solution; 

<p = 1 (6-41) 

K 3 + v K 3 - ^4 cosh r 



where we have chosen the peak position to be at r = 0. This solution is apparently meaning- 
ful only for k\ > K4, which is a valid assumption as we will see later. Next, we must solve the 
self-consistent equation (6-39) to determine the value of 9. Substituting the solution (6-41) 
into Eq. (6-40), we obtain the following equations for #3,4; 

r 3 = -^ + ^"^ arctanh^, (6-42) 



«4 Kl' V K 4 

R 4 = zM±±!i + %(5^-3. 4 ) arctanh ^_ (643) 

6K 4 K ^ ^[K 4 



In order to find solutions, we find it convenient to choose 9 and K4. Using the relation 
2^3 cos # — K4 + I, which can be derived from Eqs. (6-36) and (6-37), we express the above 
functions R3 4 as functions of the two parameters, 9 and /t 4 . And then we solve two equations 
(6-37) and (6-39) for 9 and /t 4 . (Because Eqs. (6-37) and (6-39) are invariant under the 
reflection 9 — > 71 — 9 and — > —0, we investigate only the region where 9 G [0, 7r/2].) In 
Fig. 30, we plot the numerical solutions of these two equations for n = 15; the solid line 
satisfy Eq. (6-39), and the dashed line Eq. (6-37). The solutions of (6-31) are given by the 
cross points of these lines. The number of the solutions are obtained from the following 
consideration: The parameter a varies from it/ '4 to (n + 3/4) n when 9 varies from to 
n. Since k 4 = cos 29 — sin 29 cot a, there is always one solution for every region width n 
in a. Therefore the number of solutions is either n or n + 1. All these solutions are the 
configurations that are stationary points of the perturbative functional F n [0] , which we have 
named the "asymptons" . 
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Fig. 30. The solid line satisfy Eq. (6-39), 
and the dashed line Eq. (6-37)(n=15). 

In order to see which asympton dominates the path-integral, we have evaluated the 
value of F n [(f>] at an asymptons, 4>+(9 asym £ [0, 7r/2])), and its partner ( "anti-asympton" ) 
0_ = — 0+(7r — Qasym £ [ 7r /2, 7t])- The result is the following; 



1 



n 



V im sin 9 Vi? 4 

(-)^n[0 + 



exp(— /t 4 — 2 sin 2 0)^ 



n/2 



sin a, 



(6-44) 
(6-45) 



Most notable is the factor rW 2 in F n [<p}. This agrees with the leading term of the known large 
order behavior of the coefficients, z n ~ n n ^ 2 . In this sense the existence of the asympton 
explains the non-Borel summability in the perturbation series. We note that Eq. (6-45) 
shows that for odd n the contribution of the asympton and the anti-asympton cancel. This 
is in agreement with the fact that z n = for odd n. 

The expression Eq. (6-44) shows that the leading asympton is the one with the largest 
e' K ^ 2sin2e /R 4 (= /(« 4 ,0))- Fig. 31 shows the behavior of /(k 4 ,0) in the (k 4 , 0) space. The 
solid line satisfy Eq. (6-39). On this line, the values of /(/t 4 ,0) is represented by Fig. 32. 
The maximum point is numerically 9 ~ 0.6(/t 4 ~ 0.3). The asympton, of which the value of 
is most close to that of the maximum point, is the leading asympton. 



6.2.3. Asympton in the SU(2) gauge field theory 

The analysis we have done for the quantum mechanics can be straightforwardly extended 
to the SU(2) gauge field theory. This may provide the way to predict the large order behavior 
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/(K 4 ,e) 




Fig. 31. The function f(K^,0) in the («;4,i 
space. 



e 



Fig. 32. The values of /(«4, 6) on the line that 
satisfies Eq. (6-39) 



of the model.0 The action of this model is 

S[A fl ,g] = ^Jd 4 xTr(F, u F, 
where we define the field strength as, 



[IV } 1 



[IV 



d^A u - d u A tl - ig[A^ A v ] 



(6-46) 



(6-47) 



In this model, the effective equation of motion (6-32) that determines the asymptons reads 

(6-48) 



as follows; 

lm[{d v F ia , + ig-.[Fp,,A v \) e i °\ = Q. 
We solve this equation (6-48) under the following ansatz; 



A^x) 



ST 



in* x" 

where r = logs 2 (x 2 = Y%=i x 1) anc ^ us i n g the Pauli matrices cr M , 



(6-49) 



(6-50) 



This ansatz is motivated by the fact that the instanton, anti-instanton, and instanton-anti- 
instanton pairs (valley) that shares a common center x^ = 0, can be described using this 
ansatz. The existence of interplay between perturbative and non-perturbative effects implies 



*> Actually, there is another cause of the non-summability in field theories, called renormalons. It is 
most plausible that renormalons are not isolated configurations, but a sum of contribution from the path- 
integration over a wide region of the functional space, unlike the asymptons. This is still an open question. 
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that they share similar structure. We will see that it is evident from the analysis of the 
asymptons. 

By substituting (6-49) in (6-48), we find that (6-48) reduces to the following single equa- 
tion; 

U 2 s' 



Im 



dr 



2 



- 2g 2 _s 3 + 3<M 2 - s )e itT 



0, (6-51) 



where g_ = g_/V487r 2 . It is important to stress here that this is obtained not by narrowing 
our search for the asympton down to the subspace described by the ansatz Eq. (6-49); as 
long as Eq. (6-49) is satisfied, the full asympton equation Eq. (6-48) is satisfied. Therefore 
if we find a solution to Eq. (6-49), it is a peak of the perturbative functional F n [(p] in the 
full functional space. This equation can be written as follows by using the scaled variable 

S ( T ) = l#-K r ) : 

d 2 s , dV eS (s) 



dr 2 ds 

where the effective potential V e g(s) is defined by 



0, (6-52) 



V eS {s) = ^S 2 - k 3 5 3 + ^ 4 5 4 - (6-53) 

This effective equation of s equals to that of the quantum mechanics (6-34). Therefore we 
can use the result of the quantum mechanics. This way the asymptons are identified in 
the SU(2) gauge field theory. Just as the asymptons in the quantum mechanical model were 
quite different from instantons or instanton pairs, and rather like the bounces, the asymptons 
in the SU(2) gauge field theory are not like instantons at all. It has non-zero field strength 
on a spherical shell, and look like a merged concentric instanton-anti-instanton pair. 

6.2.4. Asympton on the Valley 

From the above we have found that asymptons are quite unlike the instanton configu- 
rations in general. They are off the valley line studied extensively in the previous sections. 
On the other hand, we have seen in section 6.1 that the consideration of valley provides 
an insight into the summability properties of the perturbative series. Thus we will next 
study the behavior of the perturbative functional F n [<j)} on the valley.0 In order to simplify 
the analysis, we adopt the following approximation for valley in the double-well quantum 
mechanical model; 

s( T ) = -(1 + e- T - d/2 )(l + e T ~ d/2 ), (6-54) 
9 



*) Similar analysis were presented in Ref. p3| ). However, there was a mistake of factor two in one of the 
expression obtained in that paper. This changes the behavior of F n [4>] drastically. 
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In (6-54), if we choose d large enough, we obtain instanton-anti-instanton pair. It gives a 
reasonable approximation of the valley line parametrized by d. 

In the SU(2) gauge field theory, the valley has the tensor structure (6-49). If we substitute 
this expression into the action, we find that 



/oo 
dr 
-oo 



1-2 1 2/1 ~ \2 

2 S ~2 S {1 - 9S) _ 



(6-55) 



where the r-integration comes from the integration over the four-dimensional radius from 
zero to infinity, and g = gj\j / 487r 2 . Since the right hand side is equivalent to the double- 
well quantum mechanical model with (imaginary) time r, and the coupling constant g, the 
following analysis is also valid for the SU(2) gauge field theory by rewriting g by g. 

We first calculate the ratio of F n [A^\g n to e~ s on the valley. Using (6-27) for Eq. (6-30), 
we find it to be the following; 

Fn[<P}9 n A ■ (a cc\ 

5 — = e sin er, (6-56) 

e~ b 

where 

A = — + - - gc x + g c 2 --log— log(7msm0). (6-57) 

0C2 2 2 2q t z C2 2 

The magnitude of the perturbation coefficient is determined by e A . Therefore, we examine 

the behavior of e A for large d. The coefficients ci^ are approximated as follows for large d; 



ci = \ (d — 3 + • • 
C2 = ^[d-— + - 



2g* V 3 



(6-58) 



Using (6-58) for A in (|6-57|) , we find that 



From (6-59), we find that e A has the peak on d = 2ng 2 and that the magnitude of the 
peak is proportional to 2 n ^ 2 /y/n. We will call this peak "the valley asympton", as it is the 
asympton solution when limited to the valley line. As n increases, the valley asympton 
moves to the region of large d and grows exponentially. The numerical analysis shows that 
the main contribution of the asympton is in the following region: 

0.46ns 2 + y < d < 5.34V + y . (6-60) 

We can obtain this result by solving the equation A = numerically. The perturbative 
functional F n [<f\ on the valley is plotted in Fig. 33. 
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Fig. 33. The behavior of F n [<j>\g n /e~ s for n= Fig ' 34 The behavi ° r ° f ^=0 F n [<t>\g n &r 
30, g = 0.4. N = 30,g = 0.4. 

The numerical comparison of the partial sum, 

En=o F n[<P]g n with e~ s is plotted in Figs. 34, 
35 (for N = 30), and 36 (for iV = 50). (Note the difference in the vertical scales between 
the first two and the latter.) 



n=0 6 




0.46Ng+^L 



10 



n=0 6 




0.46Nt + ±f 10 



Fig. 35. The dashing line is e s . The solid Fig. 36. The dashing line is e~ s . The solid 
line is E^=o F n[4>]9 n for N = 30, g = 0.4. line is £^=o F n [<f>]g n for TV = 50, g = 0.4. 



From Figs. 34, 35, we find that Y,n=o F n[ ( t ) \g n deviates from e s where F n [<f>] has the large 
contribution. For d<0A6Ng 2 + ^, we find that £n =0 F n [(j>]g n (the solid line) is very close 
to e~ s (the dashed line). This means that as far as the valley is concerned the n-th order 
perturbation is valid only in the region d < 0A6Ng 2 + y. The usual perturbation suffers 
from non-Borel summability due to the valley asympton, even when limited to the valley 
line. 
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6.2.5. Discussion 

We have so far identified the causes of the Borel non-summability as asymptons and 
valley asymptons. the next problem is use this knowledge to look into the interplay of the 
perturbative series and the non-perturbative terms. In the double-well potential model, 
Aoyama and proposed a modified perturbation theory, obtained by dividing 

the functional space into two fundamental regions of the symmetry — > l/g—(j), and carrying 
out the functional integral in only one of them. The original path-integral is recovered by 
multiplying two on that integral, as the contribution from the fundamental regions are the 
same. As the order of the perturbation exceeds 0(l/g 2 ), the asympton moves out from a 
fundamental region into another. Thus the modified perturbation series is free from the 
asympton at higher orders. What remains is the anti-asympton, which have the property 
(6-45). This still renders the series, non-convergent. But due to the disappearance of the 
asympton, the cancelation for n = odd no longer holds, and the modified perturbative 
coefficient z n behaves as, 

z n ~ {-l) n n n '\ (6-61) 

for n — > oo. Thus this series is a Borel-summable series, which is quite desirable for many 
purposes. 

This method, however, works only for the finite time-interval (or finite temperature). 
When the time-interval is infinity, the boundary of the fundamental region is at infinity and 
the asympton never moves out of the original region. Furthermore, we do not know how to 
make it work in the SU(2) gauge theory. 

The valley asympton gives us a hope for possible resolution of the problem. It shows 
that (at least on the valley line) the perturbative series suffers from non-summability due 
to the configurations (valley asymptons) that live in d > 0(1) x n/g 2 region. This means 
that if we introduce a cutoff on d at 0(1), the perturbative coefficients begins to converge 
for n < 0(1)/ g 2 . Since this is the order when the usual leading term, n n l 2 g r \ is of 0(1), this 
perturbative series with the cutoff never suffers from any large terms. On the other hand, the 
region excluded by this cutoff is where the instanton and the anti-instanton is separated well. 
Thus the non-perturbative calculation (valley method) provides a quite effective evaluation 
of the path- integral. 

What is to be done is to find a way to apply this kind of the cutoff method to the whole 
region with all the virtues seen on the valley. This is a fascinating possibility, but remains 
as an open question at this moment. 
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